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COMPUTER PROGRAM FOR DESIGN ANALYSIS 


OF RADIAL-INFLOW TURBINES 
by Arthur J . Classman 
Lewis Research Center 


SUMMARY 

This report presents a computer program for the design analysis of radial-inflow 
turbines. Input design requirements are power, mass flow rate, inlet temperature and 
pressure, and rotative speed. The design variables include stator-exit angle, rotor- 
exit-tip to rotor-inlet radius ratio, rotor- exit-hub to tip radius ratio, and the magnitude 
and radial distribution of rotor-exit tangential velocity. The turbine losses, which are 
determined by an internal loss model, include those due to stator and rotor boundary 
layers, tip clearance, disk friction, and exit velocity. The program output includes 
diameters, total and static efficiencies, and all absolute and relative temperatures, 
pressures, velocities, and flow angles at stator inlet, stator exit, rotor inlet, and rotor 
exit. At the rotor exit, these values are presented at any number of radial positions up 
to a maximum of 17 . 

Presented in this report are the loss model, the analysis equations, an explanation 
of input and output, and the FORTRAN program listing and variable list. Sample cases 
are included to illustrate use of the program. 


INTRODUCTION 

The analysis of a power or propulsion system involves many repetitive calculations 
of component performance and geometry over a range of conditions . Such calculations 
are most easily and quickly done by a computer. One component of interest for small 
gas turbine systems is the radial-inflow turbine. Radial-inflow turbine geometries for 
achieving maximum static efficiency are presented as a fimction of specific speed in 
reference 1. However, there appeared to be no readily available computer program for 
performing the velocity-diagram analysis required for determining geometry and esti- 
mating performance. 



A computer program for the design analysis of radial- inflow turbines was, there- 
fore, developed. Input design requirements are power, mass flow rate, inlet tempera- 
ture and pressure, and rotative speed. The design variables include stator-exit angle, 
rotor- exit- tip to rotor-inlet radius ratio, rotor-exit-hub to tip radius ratio, and the 
magnitude and radial distribution of rotor-exit tangential velocity. The turbine losses 
include those due to stator and rotor boundary layers, tip clearance, disk friction, and 
exit velocity. The program output includes diameters, total and static efficiencies, 
and all absolute and relative temperatures, pressures, velocities, and flow angles at 
stator inlet, stator exit, rotor inlet, and rotor exit. At the rotor exit, these values can 
be presented at any number of radial positions up to a maximum of 17 . 

This radial-inflow tui'bine design analysis computer program is described herein. 
Presented in this report are the loss model, the analysis method, an explanation of in- 
put and output, and the FORTRAN program listing and variable list. Sample cases are 
included to illustrate use of the program. 


LOSS MODEL 

An important part of any turbine design problem is the estimation of losses. The 
loss model used for this analysis is a modification and extension of the model used in 
reference 1. Accounted for by this model are the three-dimensional (profile plus end 
wall) viscous losses in the stator and the rotor, the disk-friction loss on the back side 
of the rotor, the loss due to the clearance between the rotor tip and the outer casing, 
and the exit velocity loss. 


Viscous Loss 


The stator and rotor viscous losses are each expressed in terms of a kinetic- 
energy loss coefficient, which is defined as the loss in kinetic energy as a fraction of 
the ideal kinetic energy of the blade row actual flow. In terms of boundary layer param- 
eters (see ref. 2), the two-dimensional kinetic- energy loss coefficient is expressed as 


®2D 




s cos $ - ^ 


( 1 ) 


where the flow angle $ is equal to for the stator and jSg for the rotor. 

The symbols are defined in appendix A. The station designations are indicated on the 
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schematic cross section of a radial-inflow turbine shown in figure 1. An example 
velocily diagram indicating an^e designations is shown in figure 2. Substituting 

>p = Ed (2) 


and 


6 = H0 (3) 

into equation (1), introducing the surface length I, and dividing numerator and denomi- 
nator by blade- row exit spacing s yield 



(4) 


The momentum thickness per unit of surface length is assumed to be expressed as 
a function of Reynolds number as 


^tot 

I 
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(5) 


The factor C is intioduced as a convenience for modifying the loss level if desired. 

It is further assumed that the ratio of three-dimensional loss to two-dimensional loss 
is equal to the ratio of three-dimensional (blade surface plus end wall) surface area to 
two-dimensional (blade wall) surface area (see ref. 2); that is. 


®3D 



( 6 ) 


Combining equations (4) , (5) , and (6) yields 
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Equation (7) is used to determine tiie three-dimensional viscous losses in the stator and 
the rotor. Evaluation of the various terms in equation (7) is as follows. 

Flow angle . - As mentioned previously, 

= “l <«> 


and 


Energy and form factors . - The energy and form factors are obtained from equa- 
tions (7-14) and (7-13), respectively, of reference 2 using a velocity profile exponent 
of 0.2. Sufficient numbers of terms are used for each series such that additional 
terms affect the series sum by less than 0. 1 percent. The resultant equations are 


j._ \1.92 3.2 4.8 6.72/ 

1.68 2.88 4.4 6.24 

and 

H - 1-2 1.6 2.0 2.4 2,8 

1.68 2.88 4.4 6.24 


( 10 ) 


( 11 ) 


For the stator, 


4 


2 


( 12 ) 



and for the rotor, 


Q, 


_ 7 - 1/ W 
y + 1 \W, 


CT/ 


(13) 


Refe rence loss coefficients . - The reference values of momentum thickness per 
unit length and Re 3 molds number were combined into reference loss coefficients, 
which were evaluated by matching the program results with the experimental values of 
stator loss coefficient (^ 3 j) g “ 0 - 064, which was obtained from unpublished data) and 
total efficiency (t]' = 0.88) for the turbine of reference 3. This is considered to be 
representative state-of-the-art performance for a carefully designed turbine. With the 
coefficients and equal to unity, the reference loss coefficients are 

= 0.03734 (14) 

,s 



and 


tot 


\JL Re 


- 0 . 2 , 


= 0.11595 


ref, r 


(15) 


Re,ynolds number . - Blade- row chord is used as the characteristic length. For the 
stator. 


and for the rotor 


Re 


Pl^lS 


Re 


r 




jU 


(16) 


(17) 
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Expressing velocity in terms of its thix)ughflow component, multiplying both numerator 
and denominator by annulus area, and applying the continuity equation yield 


Res 


w 



cos 


(18) 


and 



(19) 


Length to spacing ratio . - For the stator, the surface length to exit spacing ratio is 
expressed as 



( 20 ) 


A stator solidity (c/s)„ is input to the program and internally adjusted to correspond to 

s 

an integral number of vanes. The stator surface length to chord ratio (l/c)^ is unity 
for an uncambered vane and obtained from equation (BIO) of appendix B for a cambered 
vane. 

For the rotor, the surface length is obtained from equation (C2) of appendix C. 

The mean spacing at the rotor exit is 

27rr« 7r(r2 v, + J 

s = ^ = 111- (21) 

n n^ 

r r 

Trailing-edge thickness to spacing ratio . - TraUing-edge thickness is specified as 
a fraction of blade height using existing turbines as a guide. For the stator, the 
trailing-edge thickness is 5 percent of the stator blade height; that is. 

The solidity c/s is an input value, while the aspect ratio h/ c is either input or com- 
puted from the stator geometry model depending on the case. 
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For the rotor, the tralling-edge thickness is 4 percent of the rotor blade exit 
height; that is, 


t = 0. 04(r2 ^ 


^2,h^ 


Dividiag equation (23) by spacing as expressed by equation (21) yields 




( 23 ) 


(24) 


Three-dimensional- to tw o-dimensional-area ratio . - Stator vane surface area 
for one passage is 





(25) 


Stator end-wall surface area for one passage is 



Expressing number of stator vanes as 

27rr, 

and substituting equation (27) into equation (26) yield 



Since 


( 26 ) 


(27) 


(28) 
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i' 


( 29 ) 


^3D 

■^2D 


A, + A 
b w 


substitution of equations (25) and (28) into equation (29) yields 


For the rotor, 



(30) 


(31) 


where A, A,„ . and A,„ , ^ are obtained from equations (C5), (C13), and (C17), 
respectively, of appendix C. 


Disk- Friction Loss 


The disk-friction loss is calculated using equations (8-7), (8-10), and (8-15) of ref- 
erence 4 with the coefficient Cjy = 0. 085. For one side of a disk, the friction loss is, 
therefrom, expressed as 





la 


(32) 


Tip Clearance Loss 

In accordance with the tip clearance discussion in reference 5, it is assumed that 
the fractional loss due to clearance is equal to the ratio of clearance to passage height 
at the rotor exit. Clearance is input to the program as a fraction of diameter. The 
loss due to clearance is, therefore, computed from 
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Ah. 


VD,av ^2,t~^2,h VtL ^2,t~ ^2,h 


/h 
= c 


2r, 




D. / To 4. - r„ 



where is the input parameter. 


Exit Velocity Loss 


The kinetic- energy loss associated with the turbine exit velocity is 


( 33 ) 


2gJ 


(34) 


ANALYSIS METHOD 


The flow analysis is one dimensional at the stator inlet, stator exit, and rotor Inlet, 
each of these calculation stations being at a constant radius . At the rotor exit, where 
there is a variation in flow-field radius, an axisymmetrLc two-dimensional analysis is 
made using constant height sectors. Simple radial equilibrium is used to establish the 
static pressure gradient at the rotor exit. 

The fluid energy corresponding to the shaft power requirement is 


^shft 


Jw 


(35) 


Disk friction and clearance losses, which are expressed by equations (32) and (33), re- 
spectively, are added to the shaft work to yield the average total fluid energy extrac- 
tion; that is. 


^VD 


= Ah’... 

av shft 


+ L + L 
df ( 


(36a) 


which can also be expressed as 
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(36b) 


Ah' 


VD, av 


^shft ^df 


1 - 


Ah 


VD, av 


The numerical values needed for evaluating are not immediately known, and an 
iteration is used starting with = 0. 

The tuibine energy transfer equation is 


Ah'™ 
VD, av 


= ^KaVu,Xa-;Z«l'^2.lV2.i 

\ i=l / 


(37) 


Substituting 


^2 i 

"2.1 = 7^’Jla 
na 


(38) 


into equation (37) and manipulating terms yield 


Ah;,^ (i 

gJ 


VD, av 


K > 

-lV 

'*' / a ‘ ’^la''u,la 


i=l 


The average change in tangential momentum is 


(39) 


<-^''u)^„=--laVla-t:E*i'^2.lV2.1 

i=l 


(40) 


As discussed in reference 5, the rotor inlet gas tangential velocity and blade speed can 
be related as 


V, 


u, la _ ^ 2 


U. 


(41) 


la 


n. 
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Substituting equations (40) and (41) into equation (39) and solving for blade speed yield 


-il/2 



The ratio is specified as input, and n^ is determined from equation (Cl 9) of 

appendix C. With known, equation (41) is used to determine Rotor inlet 

diameter is 


U. 




la 

ttN 


(43) 


Stator Exit 

The conditions at the stator exit, station 1, are computed from the following equa- 
tions. The radius ratio the stator exit angle oi^ are known input values. 

Stator loss coefficient e„ can be either input or computed from the previously pre- 
sented loss model. If the loss model is used, the stator-exit calculation is iterative 
between equations (50) to (54) . 


D 


1 



V 


u,l 



V 


1 


V 


iLi. 


sin 01 . 


(44) 


(45) 


(46) 




tan 


d. 


(47) 
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(48) 


T* =T« 


v: 


Ti = T- - 


2gJc 


y2 =-!L- 
l,id , - 

1 - 


V? . , 


vy/ (y-1) 


f \y/ (Y-i) 


Pb 


/TrV 

pi'Pi(:;r) 

\ 1/ 


Pi 


Pi = 


RT, 


h = 
s 


w 


Pl^Kld, 


V = 'tjjtl— 

er,l 


gRT^ 


(49) 


(50) 


(51) 


(52) 


(53) 

(54) 


(55) 


Stator Inlet 

Stator- inlet diameter Dq and flow angle Q!q are either input or calculated from 

the stator geometry model, appendix B. Stator height h is assumed constant, and the 

s 

following four equations are solved simultaneously for the stator-inlet, station 0, con- 
ditions Vq, Tq, Pq, and p^; 


’’^^O^sPo “o 


( 56 ) 


12 



( 57 ) 




V? 


2gJc 




P 

vy/ (y-1) 


Po = P6 - 


Po 


Pq 

RTn 


(58) 


(59) 


Rotor Inlet 

Assuming no change in total temperature or total pressure between stator exit and 
rotor inlet (i. e. , and pj^^^ = p^) , we can solve the following set of six equa- 

tions for the variables indicated on the left sides: 


'la 


u, la 


sin a 


la 


(60) 


Tla = Tia ‘ 


^a 

2gJc„ 


(61) 


Pla^Pla 


la 


\y/ (y-i) 


lay 


(62) 


^la 


_ Pla 


RT 


la 


(63) 


^rad, la ^rad, 1 


Pl^l 

Pla^la 


(64) 
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(65) 


= tan“^ — 


V 


rad, la 


At this point there is sufficient information for evaluation of disk friction loss from 
equation (32), and the calculation cycles back to equation (36b) until a convergence is 
obtained. 

The rotor-inlet conditions relative to the rotor are 


Wu,ia= Vla-^la 


( 66 ) 


^la = tan' 


-1 u, la 


V 


rad, la 


(67) 


W - ^^ad, la 
la 


cos j3 


la 


( 68 ) 


W 


T" =■: T + 
^la ^la 


la 


2gJc 


(69) 


pn \y/(y-i) 


Pla 


la' 


^la 


T 


la> 


( 70 ) 


IT' 

= V / ia 

cr, la cr, 1 j, rp, 
^la 


f71) 


Rotor Exit 

At the rotor exit, station 2, the annulus is divided into equal-height sectors, with 
the sectors being related by simple radial equilibrium. The solution at the rotor exit 
has to satisfy two input requirements, average specific work and flow. This is done by 
means of two iteration loops. The outer-loop iteration is started by setting a value for 
the mean- sector tangential velocity 2 ni' inner loop is then cycled by varying 
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the mean-sector axial velocity V « until continuity is satisfied. The average spe- 
cific work is then checked. If necessary, g m is adjusted and the iteration process 
repeated until both loops are simultaneously satisfied. 

The rotor-exit hub and tip diameters are known from the calculated rotor-inlet 
diameter and input diameter ratios. For the k sectors, therefore, the sector mean 
diameters D„ . are known. Without knowledge of the rotor-exit solution, we can then 
calculate for all i 


7TD„ .N 
TT = ^ 

c 


(72) 


-T” 

^2,1- 2gJCp 


(73) 


w 


, =w 1^ 

cr, 2,i cr, la I/ rptf 

f la 


/v; .\Y/(r-i) 

P2,id,l = Plal^' 


la, 


(74) 


( 75 ) 


The rotor loss coefficient e^ can be either input or computed from the previously 
presented loss model using mean-sector conditions. With the rotor-exit mean-sector 
tangential and axial absolute-velocity component values set by the iteration loops pre- 
viously discussed, the rotor-exit calculation proceeds as follows, starting with the 
mean sector: 


^u,2,m ^u,2,m”'^2,m 


139. ^ = 


V, 


X, 2,m 


W 


2, m 


_ ^x, 2.m 


cos /3 


2,m 


( 76 ) 

(77) 


( 78 ) 
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( 79 ) 


rp rpTf 

^2,m ^2,m 


w; 


2, m 


2gJc 


wS =_2 l“ 

2,id,m ^_- 


(80) 


W2... 


yy/ (y-i) 


P2 




(1 2,id.m 

2gJc TU 

^ p 2, my 




m 


(81) 


P2 


P2, 


m 


,m 


RT, 


2, m 


(82) 


For the other sectors, the rotor loss coefficient is assumed constant and equal to 

the mean-sector value. The angular momentum distribution (rV ) /(rV ) is 

\ u/2^i/\ ^/2,ra 

specified as input. The calculation then proceeds from the mean sector into the hub and 
from the mean sector out to the tip; thus, 


V ,. = v , 

U.2., 


/2,m 


In the following, the subscript j refers to all values of i except i = m: 


''^u,2,j ^u,2,j“^2,j 


P2, , = P2. j.l " °2, ,*l) 

^g \U2 J ^2,j±l 




p \(y-i)/y 

1 - I 

\P2,id,j/ 


(83) 


(84) 


(85) 


( 86 ) 
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( 87 ) 






w 

jSo , = sm~^ — UjAlJ. 
J w 

^2,] 


( 88 ) 




( 89 ) 


WS 


T2 , = T-^ i - 

2gJc 


(90) 


P2. 


’ ‘'^2J 


(91) 


After equations (83) to (91) are evaluated for all values of j, the mass flow rate at 
the rotor exit is calculated as 


'^calc = 


^2,t ~ ^2,h 

2k 


XV 

X ^2,i^x,2,i^2,i 


i=l 


(92) 


If ^ w, the value of 2 m adjusted, and the calculation is cycled back to 

equation (77) . When energy transfer determined from equation (37) is 

compared with the work requirement from equation (36b) and, if necessary, the value of 
Vu 2 m adjusted and the calculation cycled back to equation (76). When the flow and 
work requirements are simultaneously satisfied, the rotor-exit solution is established. 
The following parameters are then calculated: 


Ah 


VD, i 


_ 

gJ 


(U. 


la^u, la “ ^2, 


i^u. 


,i> 


(93) 


Ahi 




VD.i 


(94) 
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^cr, 2 , i ^cr, la 



V., 


otg . = tan ^ 

x,2,i 




y/ (y-1) 



^x,2.i 


cos 01 


2,i 


Performance 


There are five types of losses considered in this analysis. The disk' 
clearance losses are obtained from equations (32) and (33), respectively, 
rotor, and leaving losses are 


2gJ 


k 


L 


r, av 



i=l 


w. wf . , . - W? . 
1 2,id, 1 2,1 


w 


2gJ 


k 


L 


ex, av 



1=1 


w. vf . 
w 2gJ 


The sector ideal specific works are 


(95) 


(96) 


(97) 


(99) 


friction and 
The stator. 


( 100 ) 


( 101 ) 


(102) 
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= Vo 


Pii 


(y-l)/y 


1 - ^ 


Po> 


Vo 


1 - 


^p^i 

Po 


,(y-i)/"M 


The sector efficiencies based on velocity diagram work are 


'^VD,i 


^VD.i 

^ld,i 


Ah 


VD,i 


'VD, i 


Ah, 


id, i 


The overall velocity diagram efficiencies are 


Ah 




VD,av 


VD, av k 


z 


i=l 


w. 


Ah! 


w 


id, i 


Ah 


’^VD.av k 


VD, av 


z 


i=l 


w. 

— Ah., . 
w 


while the net overall efficiencies are 


(103) 


(104) 


(105) 


(106) 


(107) 


(108) 


^shft = ^VD, 


^^shft 


av 


Ah 


VD, av 


(109) 
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( 110 ) 


’’shft '^VD.av^, 


^shft 


VD, av 


A design parameter that is computed by the program is specific speed 


N 




(• 


J Ah! 


id, aV/ 


)3/4 


( 111 ) 


DESCRIPTION OF INPUT AND OUTPUT 

This section presents a detailed description of the program input, normal output, 
and error messages. The input and corresponding printed output for an example prob- 
lem are included for illustrative purposes . 


Input 

The program input consists of a title card and a data set in NAMELIST form for 
each case. The title, which is printed as a heading on the output listing, can be located 
anywhere in columns 1 to 78 on the title card. A title card, even if it is left blank, 
must be the first card for each case. 

The data are input in data records having the NAMELIST name INPUT. The varia- 
bles that comprise INPUT along with descriptions, units, and special remarks are pre- 
sented in the list to follow. Either SI units or U. S. customary units may be used with 
this program. Values for some of the variables in the input list are internally preset 
by the program before reading the input. Thus, if a preset value (see input list for val- 
ues) is appropriate, that particular variable does not have to be specified in the input. 

lU units indicator (preset value = 2) : 

1 - SI units 

2 - U. S. customary units 

NSTAR stator geometry indicator: 

0 - cambered vane with specified diameter ratio 

1 - cambered vane with specified aspect ratio 

2 - uncambered vane 
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ALPHAO 

ALPHAl 

CDT2 

CR 

CS 

EBARR 

EBARS 

GAM 

K 

MU 

N 

POW 

PTIN 

R 

RH2RT2 

RT2R1A 

RVIAAV 

RV2I2M(I), 
1= 1,K 

RORIA 

RIRIA 

SIGSESr 

STAR 

TTIN 

W 


stator-inlet flow angle from radial direction, deg (input only for 
NSTAR = 0 or 1) 

stator-exit flow angle from radial direction, deg 

ratio of clearance gap to blade tip diameter at rotor exit (preset 
value = 0.) 

rotor loss-coefficient multiplier (preset value =1.) 

stator loss-coefficient multiplier (preset value =1.) 

rotor loss coefficient (input only if it is desired that internal loss model 
not be used for computing e^ 

stator loss coefficient (input only if it is desired that internal loss model 
not be used for computing e„) 

O 

specific heat ratio 

number of radial sectors at rotor exit 

9 

gas viscosity, (N)(sec)/m ; lb/ (sec) (ft) 

rotative speed, rad/ sec; rev/ min 

shaft power, kW ; hp 

inlet total pressure, N/cm ; Ib/in. 

gas constant, J/ (kg) (K) ; (ft) (Ibf)/ (Ibm) (°R) 

rotor- exit-hub - to tip- radius ratio 

rotor- exit- tip - to rotor-inlet-radius ratio 

ratio of rotor- inlet angular momentum to average change in angular 
momentum (preset value = 1.) 

ratio of i^ sector to mean-sector rotor-exit angular momentum 
(preset value = 1 . for all I) 

stator-inlet- to rotor-inlet radius ratio (input only for NSTAR = 0) 

stator-exit- to rotor-inlet radius ratio 

stator solidity based on exit spacing (preset value = 1. 35) 

stator aspect ratio (input only for NSTAR = 1 or 2) 

inlet total temperature, K; °R 

mass flow rate, kg/sec; Ib/sec 
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The program inputs for two sample cases are shown in table I. Both cases are for 
the same design requirements. The first case is with SI units and the second case is 
with U. S. customary units. The first card for each case is tiie previously described 
title card. After the title card is the NAMELIST data set containing all the variable 
values. For the second case, only those variables changing in value from the first case 
need be input. The output corresponding to this sample input is presented and described 
in the following section. 


Output 

The program output consists of title headings, the input variable values, and the 
computed results. This section presents normal output. Error message output is de- 
scribed in the next section. 

Tables II and III present the output that corresponds to the sample input shown in 
table I. Table II is for the case with SI units. The top line of output is a program iden- 
tification title that is automatically printed with the first case of each data package. 

The second line is the title card message. The next three lines indicate the units, SI 
in this case, used for the different variables. 

The heading *INPUT* is followed by eight lines showing the input values used for 
this case. Identification of all the items on the output listing is self-explanatory. The 
zero value shown for stator aspect ratio indicates that this was not an input for this 
case. Values of 1.0000 are shown for the stator and rotor loss coefficients to indicate 
that the actual values are computed using the internal loss model. 

The heading * OUTPUT* is followed by the computed results. Absolute tempera- 
tures, pressures, flow angles, velocities, and velocity ratios, along with diameter, 
are shown for each calculation station. At the rotor exit, these values are shown for 
the mean diameter of each sector as well as for the hub and the tip. Additional output 
for the rotor-exit sectors include flow rate, specific work, and total and static effi- 
ciencies . Also shown are the computed loss coefficients for the stator and rotor, the 
stator height, number of stator vanes, and number of rotor blades. 

Under the heading * OVERALL PERFORMANCE* are the turbine total- to-total and 
total-to-static pressure ratios, diagram specific work, and both diagram and net total 
and static efficiencies. Also shown are the individual loss components as fractions of 
the turbine ideal work and the specific speed. 

Table HI is similar to table n except that U. S. customary units are used, as indi- 
cated by the title message as well as by the specified units themselves. Note that the 
values for all the dimensionless variables are the same in table III as in table II. 
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Error Messages 


The program contains seven output messages indicating the nonexistence of a solu- 
tion satisfying the specified input requirements. These messages are presented in this 
section, and their causes are discussed. In general, when one of these messages ap- 
pear, the program input should be checked for errors . 

(1) NO SOLUTION FOUND AFTER 100 ITERATIONS FOR CONTINUITY AT ROTOR 
EXIT - This message is caused by the program making 100 iterations through sub- 
routine CONTIN without a solution being found. There is no obvious reason for this 
situation except possibly for input error in the specification of rotor-exit angular mo- 
mentum. 

(2) ROTOR EXIT CHOKES AT InIAXIMUM MASS FLOW RATE = XXXX.XXXX - 
This message is caused by the choking mass flow rate for the rotor exit being less than 
the design flow rate specified as program input. If the input design requirements are 
correct, then possible corrective action includes reducing exit angular momentum, in- 
creasing rotor- exit- tip- to rotor- inlet- radius ratio, and decreasing rotor-exit-hub- to 
tip-radius ratio. 

(3) REQUIRED SPECIFIC WORK GREATER THAN ENERGY AVAILABLE IN GAS - 
See item (4) . 

(4) SPECIFIC WORK REQUIRED IN SECTOR XX GREATER THAN ENERGY 
AVAILABLE IN GAS - These last two messages are caused by the turbine-exit total 
temperature, average value or any sector value, being less than zero. Possible cor- 
rective action includes increasing turbine inlet temperature, decreasing turbine power, 
or increasing mass flow rate. 

(5) REQUIRED STATOR IDEAL KINETIC ENERGY GREATER THAN ENERGY 
AVAILABLE IN GAS - See item (6) . 

(6) ROTOR IDEAL RELATIVE KINETIC ENERGY REQUIRED IN SECTOR XX 
GREATER THAN ENERGY AVAILABLE IN GAS - These last two messages are caused 
by the computed ideal energy required by the stator or the rotor being greater than that 
available from an infinite expansion of the gas . The probable reason for this condition 
is an error in the computed or input loss coefficient. 

(7) THE PROGRAM CAN NOT FIND A SOLUTION SIMULTANEOUSLY SATISFYING 
CONTINUITY, RADIAL EQ. , AND THE LOSS MODEL AT THE ROTOR EXIT - This 
message is caused by rotor- exit relative velocity, as determined by radial equilibrium 
and radial distribution of loss, being less than its tangential component, as determined 
primarily by blade speed. Corrective action includes decreasing rotor-exit tip radius 
and increasing rotor- exit-hub to tip radius ratio. 
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PROGRAM DESCRIPTION 


This computer program consists of main program RIFTUD (^dial hiflow Turbine 
Design), blade loss coefficient subprograms EFFIC, SIMPSl, SHUB, and SHUB2, and 
rotoi>exit ccmtinuity subprograms CONTIN and PABC. The entire program is written 
in FORTRAN IV language and has been run on both an IBM 7094 and a UNIVAC 1110. 
Running time on the UNIVAC 1110 is about 1 second per design case for a five-sector 
design, hi this section, the functions of the main and subprograms are described, the 
program variables are defined,, and the program listing is presented. 


Main Program RIFTUD 


Main program RIFTUD performs all input and output operations as well as all of 
the flow analysis as presented in the section ANALYSIS METHOD. 

Program variables . - The variables used in RIFTUD are defined in terms of the 
following s 5 rmbols, which are defined in appendix A; 


ALPHAA 

Ofa (degrees) 

BETA2(I) 

^2 1 (degrees) 

ALPHAO 

qJq (degrees) 

BETIA 

(radians) 

ALPHAl 

(degrees) 

BET2(I) 

^2 1 (radians) 

ALPHA2(I) 

Q!„ . (degrees) 

BET2M 


ALPHO 

Oq (radians) 

CDT2 


ALPHl 

01 ^ (radians) 

CHRD 

C 

ALPHIA 

01^^ (radians) 

CN 

Cn/tt 

ALPH2(I) 

Qf„ . (radians) 

CP 

"P 

ALPO 

Q!q (degrees) 

CPOW 

Cp 

ALPIA 

previous value of 

CR 

Cr 

ALSTG 

<P 

cs 


ALUNC 

cam 

DELVX 

increment in 

BETAIA 

iSfa (degrees) 
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DHIDS(I) 

^^d, i 

ENT 

DHIDSA 

^^id, av 

ER 

DHIDT(I) 

^Hd, i 

ES 

DHIDTA 

^^Id, av 

ESI 

DHSHFT 

^^shft 

ETAS 

DHVD(I) 

'^VD,i 

ETASV 

DHVDAV 

'^^VD, av 

ETASVD(]) 

DOR 

ISO/tt 


DVUAVC 

(dv ) 

ETAT 


^''2,av, calc 

ETA TV 

DVUTOT 

w(dV ) 

' ^^2,av, calc 

ETATVD(I) 

DVU2AV 

(dv ) 



' ^2, av, required 

EX 

DO 

Do 

F2(I) 

D1 

Di 

G 

DIA 

Dla 

GAM 



HR 

D2(I) 

^2,i 

HS 

D2M 

Do 

2, m 


EBARR 


I 

e^ input value 


EBARS 

e_ input value 
s 

ID 

EL 

2/k 

II 

EN 

n^ 



r 

IND 

ENS 

“s 


ENST 

truncated integer value of 

1ST 


Ug before roundoff 


trancated integer value of 
Hj, before roundoff 

®r 

previous value of e^ 

^shft 
’^VD,av 
^VD, i 


’Ishft 

^VD,av 

^VD,i 
y/ (y - 1) 
'"2,i 

g 

y 


dummy index, usually re- 
ferring to sector number 

sector number increment 

sector number index for 
previous sector 

solution indicator for 
rotor-exit continuity 
calculation 

iteration counter for 
stator calculation 


25 


ISTT 

iteration counter for 

PSO 

Po 


stator calculation 

PSl 

Pi 

ITER 

iteration counter for 



rotor-exit calculation 

PSIA 

Pla 

lU 

units indicator - see sec- 

PS2(I) 

P2,i 


tion Input 



J 

J 

PTIN 

Po 

K 

number of sectors 

PTRIA 

Pla 

KK 

K +2 

PTR2(1) 

P’^i 

KPl 

K + 1 

PTR2ID(I) 

P2, id, i 

K1 

indicator for print 



control 

PTO 

Po 

LC 

L 

c 

PTl 

Pi 

LCDH 

Lc^^VD, av 

PTIA 

Pla 

LCDHIS 

Lc/^id,av 

PT2(I) 

P2.i 

LLDHIS 

^ex, av^'^id, av 

Q 

(Wp)2^av 

LRDHIS 

^r, av^'^id, av 

R 

R 

LSDHIS 

^s^'^id, av 

RHOTO 

P’o 

LW 

^df 

RHOX 

previous value of 

LWDHIS 

Ld/Ahid^ av 

RHOO 

Po 

LWX 

previous value of 

RHOl 

Pi 

M 

mean sector index 

RHOIA 

Pla 

MU 


RH02(I) 

P2,i 

N 

N 

RH2RT2 


NSTAR 

stator geometry 




indicator - see section 

RLOSS 

^r, av 

PI 

Input 

RSTG 

^cm 

IT 


POW 

P 

RT2R1A 

^2./’^la 
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RVIAAV 


TTR2<:0 




TTO 


RV2IM(I) 

TTl 

0 

T* 

RV2I2M(I) 


TTIA 

Tia 

RORIA 


TT2(^ 

^2,i 

RIRIA 


TT2AV 

T^av 


k 

T/ \ // \ 1 / 

UlA 

Ula 

SIGRV 

EH.i/K.J/’' 

i=l 

UlASQ 


SIGS 

(c/s)g 

U2(D 

^2,i 

SIGSIN 

input value of (c/s) 

s 

VCRl 


ss 

s 

VCRIA 

^cr, la 

SSPD 

^sp 

VCR2(I) 

^cr,2, i 

STAR 

(h/Og 

VOVCRA 


TGJCP 

2gJc_ 





VOVCRO 

TITLEXI) 

input/ output array for 



title card message 

VOVCRl 


TSPR 

Po/P2,i 

VOVCR2(l) 

(V/Ve4_^ 

TSO 

To 



VRO 

^rad, 0 

TSl 

Tl 





VRl 

^rad, 1 

TSIA 

Tia 

VRIA 

^rad, la 

TS2(^ 

'^2,i 


TTIN 

T|) 

VUOUIA 

(V«),a 

TTPR 

Po/P^i 

VUO 

Vo 

TTRIA 


VUl 

Vi 
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VUIA 

Via 

WCR2(I) 

"'cr.2,1 

VU20 

previous value of g m 

WGIV 

input value of w 

VU2(I) 

V.,. 

WOWCRA 


VXVCR 


WOWCRM 


VXVCRP 

previous value of 

WOWCR2(I) 

VX2(I) 

Vv 9 i 



VO 

VI 

VIA 

X, ^ , 1 

Vo 
Vl 
. ^la 

WUIA 

WU2(I) 

WlA 

V la 

V,2,i 

Wla 


W2(I) 

W2,i 

VIIDSQ 

V2(I) 

V^ . , 
'^1, id 

^2,i 

W2IDSQ(I) 

<id,i 

V2LOSS 

W2TOT 

"'calc 

ex, av 

Z1 

<T/T')i,,d 

W 

w 

Z2 

<T/T")2,id,m 

WCALC 

^calc 

Z3 

(T/T")2.1d,J 

WCRIA 

Vr,la 



Program listing . - The FORTRAN listing for main program RIFTUD is as follows: 


REAL J,N,MU,LCDH,LW,LWX,LC,LCDHIS,LWDHIS ,LLDHIS ,LSOHIS,LRDHIS 
DIHENSION D2<17),U2U7»,TTR2«17>,WCR2I17),PTR2IDI17»,VU2 J17), 
lWU2U7),DHVDI17»,TT2J17),VCR2<17JtBET2ll7»,W2ll7),TS2(17),PS2a7)» 
2U2IOS0ll7),RH02J17),F2C17»,RW2IMa7J,VX2«17),ALPH2«17)»V2<17», 
3PT21 171 ,PTR2 U 7 > , VOVCR2 1 1 7 > , ttOWC R2 1 17 ),DHIDT « 1 7 ) ,DHIDS <1 7 ) , 
4ETATVDU7I ,E T ASVD 1 1 7 I , AL PhA2 1 1 7 ) ,BE TA2 U7 ) 

DIMENSION RV2I2MI17»,TITLE«13) 

COMMON/E FF/6AM,V0WCR1»W, PI, STAR, MU,D1,ALPH1,NSTAR,ALPH0, ALUNC,SIGS 
1 ,00,HS,CS,W0WCRM,01A,02M,HR,BET2M,EN,RH2RT2,CR 
NAMELIST/INPUT/R,GAM,ALPHAO, ALPHA1,N,TTIN,PTIN,W,P0W,RV1 AAV,R1R1 A, 
IRORl A,RT2R1A,RH2RT2,C0T2,MU,RV2I2M,K ,IU,EBARS,EBARR,STAR,NSTAR 
2,CS,CR,S1GSIN 
EBARS=1 .0 
E6ARR=1 .0 
CS = 1 .0 
CR=1 .0 

SIGSIN= 1.35 
ES=0.0 
RV1AAV=1.0 
CDT2=0,0 
DO 1 1=1,17 
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o n 


1 RV2I2H(I}rl.O 
IU=2 

ITERrO 

PI = 3. 14159 

D0R=57.2958 

RORlArO.o 

ALPHAD=999,99 

WRITE! 6t100> 

100 F0RWAT!1H1,39X,54HRADIAL INFLOW TURBINE VELOCITY DIAGRAM DESIGN AN 
lALYSISI 

2 READ !5»101> TITLE 

101 F0RMATI13A6) 

IFIITER.GT.O) WRITE 16,1011 ) 

1011 FORMATCIHI) 

WRITEf6,102l TITLE 

102 FORMATCIH ,13A6> 

READ 15, INPUT) 

IF INSTAR,EO.D» STAR=0.0 
IST = 0 

VXVCRPrO.O 

DELVXzO.Ol 

ALPH03ALPHA0/D0R 

ALPHlrALPMAl/DOR 

IFCNSTAR.NE.2) AL S T 6 - C AL PhQ+ AL PH 1 ) / 2 • 
EN=PI/30.*C11D.-ALPHA1)*TAN(ALPH1) 

ENTtAINTIEN) 

IF CEN-ENT.LT..5 ) EN-ENT 
IFIEN-ENT.GE . ,5) ENrENT^l, 

VU0UlAri.-2./EN 
WGI V=W 
KPlrK^l 
WRITE <6,103» 

103 FORMAT I 39H0THIS OUTPUT IS IN THE FOLLOWING UMTS./ 

1 129H0TEMPERATURE PRESSURE GAS CONST ROT SPEED MASS 

2FL0W POWER VISCOSITY VELOCITY SPEC rtORK DIAMETER 

3 ANGLE! 

GO TO C3,4),IU 

3 CPOWriOOD. 

J-l . 

6ri , 

CNZ2 , 

PT0=PTIN*10000. 

WRITE <6,104! 

104 FORMAT! 3X6HKEL VI N6X7HN/ SO C M 4 X8 H JL S / KG-K 5X 7H R A D /SE C 5 X 6HK G /SE C8 X 2 HK 
1W6X1DHN-SEC/S0 M 5X 5HM/SE C6X 6 H JL S / G M8 X2HC M8 X 7HDE GRE E S ! 

60 TO 5 

4 CP0W-55C. 

Jr777.649 
G=32,174 
CN=60./PI 
PTQ=PTIN*144 • 

WRITE <6,105! 

105 F0RMAT(12H DEG R A NK I NE 2X 8HL B / SO I N 4X 8HB T U/L B -R 5 X 7HRE V /Ml N 5X 6 HL B /SE 
1C8X2HHP7XPHLB/FT-SEC4X6HFT/SEC6X6HBTU/LB8X2HIN8X7H0E6REES! 

5 DHSHFTrCPOW/J*POW/W 
EX=GAM/!6AM-1.! 

CP=EX*R/J 

TGJCP=2.*6*J*CP 

WRITE INPUT VALUES 
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non 


WRITE ( 6*106 >TT IN, ALPHAO *EB ARS * R * PTI N* ALPHA 1 ,E8 A RR *6 AM ,N , STAR * 
1RV1AAV,MU*W*K,POW,ROR1A 

106 F0RHAT(8H0*INPUT*/13H INLET TEMP = ,F10.4 ,7X , 18HSTAT0R IN ANGLE =, 
1F6.2,6X,19HSTAT0R KE LOS COEF=, F6 . 4 , 5X , 1 4HG A S CONSTANT =F8.4/13H 1 
2NLET PRESS=,F10.4,7X,18HSTAT0R EX ANGLE =,F6.2 ,6X , 19HR0 TOR KE LOS 
3S C0EF=,F6.4,5X,14HSPEC HI RATI0=*F6.4/1 3H ROTAT SPEE0=,F10.3,7X,1 
48HSTAT0R ASPECT RAT=,F6.4, 

4 6X,19HR0T0R IN/OEL RVU =*F6.4 *5X ,14HVISC0SIT 

5Y =,E10.4/13H MASS FLOW = * F 1 0 . 4 * 7X * 1 IHOI AM R AT lOS , 19X * 1 9HR0T0R 
6 EX RAO SECTS=,I2/13H SHAFT PONE R= ,F 10 . 3 ,8X , 17HST A T IN/ROT IN =,F 
76.4 ,6X,19HR0T EX SECT/MN RVU=1 
K1=K 

IF(K.GT.IO) KlrlO 

WRITE (6*107) (RV2I2M(I)*l=l*Kl> 

107 F0RMAT(1H^*78X*F4.2*9F5.2) 

WRITE (6*108) RIRIA 

108 F0RMAT(31X*17HSTAT EX/ROT IN r,F6.4) 

IF(K.GT.IO) WRITE (6*109) ( R V2 12M ( 1 ) , 1 = 1 1 * K ) 

109 FORMAT) 1H^,77X*5F5. 2) 

WRITE (6*110)RT2R1A*RH2RT2*CDT2 

110 FORMAT ( 31X* 17HR0T EX TP/ROT IN = * F6 . 4 /31 X * 1 7HR0T EX HUB/TIP =F6.4* 
16X*19HCL HT/ROT EX TIP 0=*F6.4) 

LC0H=C0T2*2./(1 .-RH2RT2) 

LW=3.0 

6 0HVDAV=(0HSHFT^LW)/(1.-LC0H) 

ISTTrn 

TT2AV=TTIN-0HV0AV/CP 
IF( TT2AW.lt. 0.0) GO TO 203 
U1 AS0=6* J/VU0U1A4DHV0A V4RW1A AV 
U1A:S0RT(U1AS0) 

VU1A=U1 A“»VUOUl A 
D1A=CN*U1 A/N 
C 

C STATION 1 - STATOR EXIT 

C 

01=R1R1 A*D1A 
VU1=VU1A/R1R1A 
V1 = VU1/S1N( ALPHI ) 

VR1=VU1/TAN( ALPHl ) 

TTlrTTIN 

TS1=TT1-V1*W1/TGJCP 
25 ES1=ES 
ES=E8ARS 

IF(EBARS.CO. 1 .0. AND. IST.EC.O ) ES = .05 
IF(EBARS.EO.l.O.ANO.IST.GT.O) CALL EFFIC (ES*ER*1) 

IST=IST*1 

ISTT=ISTT*1 

7 VlIDSO=Vl*Vl/( l.-ES) 

21 = 1 .-WHOSO/ TGJCP/TTl 
IF(21.LE.0.0) 60 TO 204 
PS1=PT0*Z1**EX 
PTirPSl* (TTl/TSl )**EX 
RH01=PS1/R/TS1 
HS=W/RH01/WR1/P1/01 

WCR1=S0RT(2.*GAM/(GAM»1. )«G*R*TT1 ) 

W0WCR1=V1/WCR1 

STATION 0 - STATOR INLET 
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IMNSTAR»GT.0I go to 71 
DOr01A*RORlA 

CHRDrSORT fDO^+2^Dl ♦♦2-SORT C I D0*^2 ♦ 0 1 ♦♦a I ♦♦ 2- ( D0^^2-D 1 ^^2 ) ♦♦ 2 / 
ICOSf ALST6)^^2) )/2* 

STARrHS/CHRO 

ALUNC = ALPH1-ALPH0-A C0S« ID0^^2+D1^^2-4.^CHRD^^2)/2./DD/D 1 ) 

7D CONTINUE 

IFIIST.EO.I) SIGS^SIGSIN 
SSrCHRD/SIGS 
ENSr PI^Dl/SS 
ENSTrAINTIENS) 

IFCENS-ENST.lt. . 5) ENSrENST 
IFCENS-ENST.GE..51 ENSrENST+1. 

SSr PI+Dl/ENS 
SlGSrCHRD/SS 

IFCEBARS.NE.l.O 60 TO 26 
IFC ISTT.EO.l ) GO TO 25 
IFCABSIES-ES1)/ES.GT..0001> GO TO 25 
26 TTO^TTIN 

RHOTOrpTQ/R/TTO 

RHOOrRHOTO 

9 VO=W/PI/DO/HS/COS(ALPHOl/RHOO 
TS0=TT0-V0^^2/TGJCP 
PSO-PTO^CTSO/TTOI^^EX 
RHOXrPHOO 
RHOOrPSO/R/TSO 

IFCABSCRHOX-RHOO)/RHO0.GT..00Ol ) GO TO 9 
VU0=:V0^S1NCALPHQ> 

VROrVO^COSCALPHO) 

VOVCROrVD/VCRl 

C 

C STATION lA - ROTOR INLET 

C 

ALPHlAr ALPHl 
8 VRl ArVUl A/TANCALPHIA ) 

VlArvUlA/SINCALPHlA) 

TTlArTTl 

TS1AZTT1A-V1A*V1A/TGJCP 

PSlArPTl^CTSlA/TTn^^EX 

RHOl AZPSIA/R/TSIA 

VRl AzVRl^RHOl/RHOl A^RlRl A 

ALPl AZALPHIA 

ALPHlAr A TAN C VUIA/VRIA ) 

IFIABSf ALPHlA-ALPlA).GT..ul/DOR) 60 TO 8 
LWXZLW 

LMZ.0Q61 /G/J^RHOlA^Ul A**3^01 A^^2/W/ C RHOlA + Ul A *d lA/MU)+^.2 

IFC ABSfLU-LUX)/LW.GT..OOQll GO TO 6 

PTl ArPTl 

WUIAZVUI A-UI A 

BETIAZATANCWUIA/VRIA) 

WlArvRl A/COSCBETIA ) 

TTRl AZTSI A*W1A^^2/T6JCP 
PTR1A=PS1A^CTTR1A/TS1A)^^EX 
VCRl AZVCRl 

WCRlArvCRlA*SORTCTTRlA/TTlA) 

VOVCRAZVIA/VCRI A 
WOUCRArwiA/UCPlA 
C 

C STATION 2 - ROTOR EXIT 

C 
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M=<K*3»/2 

KK=K*2 

02<KK»=D1«*RT2R1A 
02U »r02<KK)*RH2RT2 
HR=ID2IKK)-D2U> )/2. 

EL=HR/FLO*T<K) 

SIGRV=0.0 
00 12 1=1, KK 

IFd.EO.l.OR.I.EO.KKI 60 TO 11 
D2< 1 )=02m*<2.*FL0AT«1 1-3. >*EL 
RV2IMU)=RV2I2Md-l» 

SI6RW=SI6RV-»RV2IM<I l/FLOAT IK ) 

11 U2ll»=02m*N/CN 

TTR2II )=TTR1 A4IU2II )**2-UlAS0»/TGJCP 
WCR2III=WCR1A*S0RTITTR2I IJ/TTRIA ) 

12 PTR2IDI I )=PTRlA*ITTR2II)/TTRlA»**rX 
RV2IMI1 J:RV2IMI2)-|RV2IHI3»-RV2IMI2» >/2. 

RV2IMIKK )=RV2IMIK^1)-IRV2IMIK )-RV2IMIK^l ) 1/2. 

0VU2AV:01A*VU1A*I1.-1./RV1AAV» 

VU2IM»=0VU2AW/02IH1*RV2IMIM)/SIGRV 

ITER=1 

13 WU2 (M)=VU2IH )-U2 IHI 
0HW0IM)r|UlA*WUlA-U2IM»*VU2(M) J/G/J 
TT2IH)=TT0-DHV01H)/CP 
IFITT2IH»,LT,0.0) GO TO 2£i5 
VCR2(M»rvCRlA*S0RTITT2IM)/TTlA> 

IFlITER.EO.n VXVCR = 0.5 

IN0 = 1 

1«* VX2IM)=VXVCR*VCR2IH) 

BET2IM»=ATANIUU2IM)/WX2CM> ) 

U2IM>rwX2IM»/C0SIBET2IM>) 

TS2IM»=TTR2IM)-W2(M)**2/T6JCP 

D2H=02IM) 

BET2MZBET2IM) 

UOUCRMr K2IM)/WCR2IM» 

ER=EBARR 

IFIEBARR.EO.l .0) CALL EFFiC IES,ER,2> 

15 U2IDS0(M)=W2««)**2/U.-ER) 

22=1 .-W2I0S0IM)/TGJCP/TTR2IMJ 
IFI 22.LE .0.0 » GO TO 206 
PS2IM)rPTR2ID(M)*Z2**EX 
RH02M>=PS2IM)/R/TS2IM» 

F2(M)=P1*02IM)’*EL*RH02IH)*VX2IM> 

IFIITER.EO.il UCALC=F2IM)*FL0AT(K) 
lFIITER.EO.il GO TO 51 

6 7 I=H 
ID = -1 

16 11=1 
I=I*1D 

VU2 III=VU2IMI*RV2IHIII*D2(MI/D2III 
UU2I1I=VU2(1I-U2I1I 

PS2 1 1 l=PS2III l*RH02 III l/G* I VU2 I II I **2/02 III l•^VU2^ I I **2/0 2 1111/2. 
1021 I I-02I1I I I 

W210S0 II |rT6JCP*TTR2II I* II.- IPS 2 II I/PTR2IDI I ll**ll./EXII 

23 = 1 .-W2IDS0 II I/TGJCP/TTR2II I 

1FI23.LE.0.0I GO TO 207 

W2III =SORTIW2IOSO 111*11. -ERII 

IFIW2I1I.LE.ABSIWU2IIIII GO TO 66 

BET2III=A SINI WU2I 1 1 /W2I 1 1 I 

VX2III:W2III*C0SIBET2IIII 


n n n n n o 


T5ZIII=TTR2II )-W2« II^*2/T6JCP 
RH02U irPS2II )/R/TS2II > 
iriI.EO.n GO TO 17 
IFCl.eo.KK) GO TO 18 
F2II >rPI*02f I )*EL*RH02f I )*VX 2 II) 

60 TO 16 

17 I=H 
ID=1 

60 TO 16 

18 W2T0Tr0,0 

00 19 I=2»KP1 

19 «2T0TrW2T0T^F2II > 

WCALC=:tf2T0T 
IFfK.GT.l) 60 TO 51 

69 1FI0VU2AV .EO. 0*0160 TO 21 
DVUTOTzO.O 
DO 20 1=2, KPl 

20 OVUTOT = DVUTOT*F2(I )^VU2U ) 

0VUAVC=DVUT0T/W 

IFI»BSCOVUAVC-DVU2AV)/DVU2AV*LE..OO01) 60 TO 21 
VU20=VU2!M> 

VU2(MlrVU20*DVU2AV/DVUAVC 
GO TO 13 

21 V2L0SS=D.D 
DHIOTA=0.0 
OHIDSA=0.0 
0 = 0.0 
RL0$S=0.0 

00 23 1 = 1, KK 

OHVOII > = «U1 A^VUl A-U2 « I )*VU2 ( I ) ) /G/ J 
TT2CI)=TT0-DHVD« D/CP 
IFI TT2U ).LT.0.0) GO TO 208 
VCR2CD=VCR1A*S0RT<TT2CD/TTXA) 

ALPH2f 1 l=ATAN( VU2f I >/VX2 a ) ) 

V2I D = VX2tD /COSI ALPH2I I > ) 
PT2CI)=PS2(I)*ITT2<I)/ TS2II))**EX 
PTR2 ID=PS2f I D(TTR2«D/ TS2(D)*^EX 
V0VCR2d )=V2«I )/VCR2II 1 
W0UCR2( 1 >=U2 (1 >/WCR2 (I ) 

SECTOR PERFORMANCE 

DHIDTI I l=CP*TTO*« l.-IPT2(D /PTO I ♦♦ <1 . /EX ) > 
OHIOSII >=CP*TTO*l 1 .-IPS2U)/PT0)** C 1 ./EX ) > 

IFC I.EO.l.OR.I.EO.KK) GO TO 22 

V2L0SS = V2L0SS*F2(D/W*V2« D**2/2./6/J 

RL0SS = RL0SS«F2(1 )/U*fU2I0S0(D-U2C D^*2)/?./G/J 

0=0^F2II I/RH02CI ) 

0HlDTA=0HIDTA>F2CII/y*0HIDTCI ) 
0HlDSA=0Hl0SA^F2(Ii/W*0HIDSf II 

22 ETATVOII )=OHVOCI I/OHIDTIII 

23 ETASVOII)=DHVOID/OHIDS(D 

OVERALL PERFORMANCE 
ETATV=OHVOAV/OHIDTA 

etasv=ohvoav/dhiosa 

ETAT=0MSHFT/0HI0TA 
ETAS=0HSHFT/0HI0SA 
SSPO=N^SORTIO l/f J^DHIOTA |♦♦•75 
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non 


LC=LCOH*DHWDA V 
LC0H1S=LC/DHIDSA 
LWDHIS=L«/DHI0SA 
LL0H1S=V2L0SS/0HI0SA 

LS0H1S=« VUDS0-V1*W1 »/2,/fa/J/DHIDSA 
L<?OHIS = RLOSS/DHICSA 
GO TO I30,32»,IU 
3J PT0=PTO/10O00. 

PS0=PS0/10000. 

PTlrPT 1 /lOOOa. 

PSI=PS1/10000. 

PTl ArPTlA/lOCOD. 

PSlArPSlA/lOuOO. 

PTRl A=PTR1A/1C0DC. 

oo=oc*ioc. 

01 = 01 * 100 . 

01 A=01A*1C0. 

HS=HS*100. 

DO 31 1=1, KK 
D2« 11=02(11*100. 

PT2I I )=PT2(I 1/10000. 

PS2(11=PS2(I 1/10000. 

PTR2 (I 1 = PTR2 (I 1/lOCDO. 

31 DHVOd 1=DHV0(I 1/1000. 

OHVOAV=DHVDA V/1000. 

GO TO 36 
?2 PT0 = PT0/l*t«». 

PSQ = PS0/1*I4. 

PT1=PT1/144. 

PS1=PS1/144. 

PTl A=PT1A/144. 

PSl A=PS1A/144. 

PTR1A=PTR1A/144. 

00=00*12. 

01 = 01 * 12 . 

D1A=01A*12. 

HS=HS*12. 

DO 33 1 = 1, KK 
02(11=02(11*12. 

PT2(I1=PT2(I 1/144. 

PS2 ( I 1=PS2( I 1/144. 

33 PTR2(I1=PTR2(I 1/144. 

URITE calculated VALUES 
36 WRITE (6,1111 

111 F0RMAT(9H *OUTPUT*/ 26 X , 2 ( 3H ABS , 6X 1 , 1 7X , 3 ( 3H AfiS , 5X 1 , 9X ,4 HREL ,4(5X 
1 , 3HREL 1/17X,4HDI A-,2(4X,5 hT0TAL 1,4X,2(6HSTATIC,3X1 , 4HF LO W , 3X , 5HVEL 
20-,4X,4HCRIT,3X,5HBLA0E,2(4X,5HT0TAL 1 , 4 X , 4HFL0W , 3X , 5HVEL0- , 4 X , 4HCR 
3IT/16X,5HMETER,2 (5X,4HTEMP,4X,5HPRESS1 ,3X,5HANGLE,4X ,4HCITV ,9H VE 
4L RAT,7H SPE ED , 5 X , 4 HTE HP , 4 X , 5HPRE SS , 8H ANGLE ,4X , 1 3HCI TY VEL RA 

571 

ALP0=ALPH0*D0R 

ALPHAA=ALPH1 A*D0R 

BETA1A=BET1A*D0R 

DO 34 1=1, KK 

ALPHA2( I 1=ALPH2( I1*00R 

34 BETA2(I1=BET2(I1*D0R 

WRITE (6,112 100, TT0,PT0,TS0,PS0, ALP 0,V0,VOVCR0 
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112 F0RNATI13H STATOR INLET , 2 < f 9 , 3 , F9 , 2 » , F9 . 3 , F7 . 2 , F9 . 2 , F7 . 3 ) 
WRITE16*1121 > EStHStENS 

1121 FORHATOIH LOSS COEFr ,F6 .4 , 6 7X , 1 IHST ATOR HGT = , F 7 .«* , 1 7H , NU MBE R OF v 
1ANES=*F5.1) 

WRITE (6,113I01,TT1,PT1,TS1,PS1»ALPHA1,V1,V0VCR1 

113 FORMATI13H STATOR EXIT , 2 i F 9. 3 , F9 , 2 I , F9 , 3 , F7 .2 , F9 . 2 ,F7 .3 » 

WRITE (6*114)DlA,TTlA,PTlA,TSlAtPSlA,ALPHAA,VlA,V0VCRA,UlA,TTRlA, 
1PTR1A,BETA1A ,W1A,W0WCRA 

114 FORMAT! 13H0R0T0R INLET , 2 ! F 9 . 3 ,F9 . 2 ) F9 . 3 ,F 7 . 2 ,F9 .2 , F7 . 3 , 2 F9 . 2 , 
lF9.3*F8.2»F9.2,F7.3i 

WRITE! 6*1141 1 ER 
1141 F0RNAT!11H LOSS C0EF=,F6.41 

WRITE !6*115l!02!Il*TT2!I)»PT2!II*TS2(II*PS2m ,ALPHA2II» ,V2!I > , 
1V0VCR2!I1*U2!I I *TTR2!I I *PTR2!I> «BETA2(I ) *W2 1 1 ) .W0WCR2 (I) , 1 = 1 *KK> 

115 FORHAT!13H ROTOR EXIT . 2 !F9 , 3 ,F9 . 2 I F9 . 3 ,F 7 . 2 *F 9 . 2 , F 7 . 3, 2F 9 . 2 , 
1F9.3,F8.2,F9,2,F7.3/!13X2!F9.3,F9.2 » F9 , 3 ,F 7 . 2 ,F 9 .2 , F7 . 3 , 2 F9 . 2 * 
2F9.3«F8.2*F9.2,F7.31 ) 

F2!l>=0.0 
F2!KKI=a.O 
WRITE!6,1161 ) 

1161 F0RMAT!1H ) 

WRITE!6*116I EN 

116 F0RHAT!26X,4HMASS,5X,4HD1AG/17X ,4HDIA-,5X ,4HFL0W,5X ,4HsPEC, 4X ,5HT0 
lTAL*3X,6HSTATlC/16X,5HMETtR,5X,4HRATE,5X,4HU0RK,2C4X,5HEFFlC» ,3X, 
223HMUHBER OF ROTOR BL A OE S= * F 5 . 1 ) 

WRITE i6*117>iD2II),F2!I ),OHVO!I »,ETATVD(I »,ETASVD«I ) ,1=1 ,KK ) 

117 FORMAT! 13X,F9,3,F9.4 ,F9,3,2F8.3 I 
TTPR=PT0/PT2IM» 

TSPR=PT0/PS2!M» 

WRITE !6,118>TTPR,LS0HIS,TSPR,LRDMIS,DHVDAV,LWDHIS,ET ATV ,LCDHIS, 

ietasv*llohis*etat*etas*sspo 

118 FORMAT !22HO*OVERALL PERF ORM A NCE*/48 X , 20HL0 SS/ I DE AL T-S DEL H/6X,24 

IHTOT-TOT PRESSURE RATIO = ,F8 ,4 , 1 3X , 1 OHS T AT OR = , F 6 ,4 /6X , 24 HTOT -ST 

2AT PRESSURE R A T I 0= ,F8 . 4 , 13X , lOHRO TOR = , F 6 , 4 /6X , 24HDI A G AV6 SPEC 

3IFIC WORK =*F8.4,13X*10HW1NDAGE = * F6 .4 /6X ,2 4H0 I AG TOTAL EFFICIENC 

4T =,F8.4,13X,10HCLEARANCE=,F6,4/6X, 24HDIAG STATIC EFFICIENCY =,F8 
5*4,13X,10HEXIT KE = , F6. 4/6X , 24HNE T TOTAL EFFICIENCY =,F8.4/6X, 
624HNET STATIC EFFICIENCY = , F8 . 4 , 8 X , 1 5HSPE C I F IC SPEED=,F7.3) 

60 TO 2 

51 lF!INO.GE.6.ANO.ABSIW6IV-UCALC)/W6IV.LE..OOOn GO TO 65 
CALL CONTIN! VXVCR*WCALC,IN0,1*W6IV* .05) 

IF! INO-10)14, 61*61 

61 IF!ITER-1)62,62,63 

62 WXVCR=.9 
IF!K.E0.1) go to 63 
IN0 = 1 

ITER=ITER*1 
60 TO 14 

63 IF!IN0-10)201, 201*202 

65 IF! ITER-1)66*66*69 

66 ITER=ITER*1 
IN0 = 1 

60 TO 67 

68 1F!VXVCR.6E.VXVCRP) 60 TO 72 
VXVCR=VXVCRP-OELVX 
DELVX=0ELVX/10. 

IF!OELVX.LE. .00001) GO TO 209 

72 VXVCR=WXWCR^DELVX 
VXVCRP=VXVCR 
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71 CHRO=HS/STAR 

iriNSTAR.tO.2) 60 TO 76 

RST6 = tCHR0#C0SI ALSTGI*S0RT<CHR0**2* fCOSI ALST6I**2-1. »/2. 
00=1 S0RT<CHR0**2/2.*2.*RST6**2-01**2/«*, ) J*2. 

ALUNC=ALPH1-ALPH0-A COSJ (00**2*01**2-«I**CHRO**2I/2, /DC/01 ) 

60 TO 70 

76 ALPHO=ATAN< SIN C A LPHl » / 1 CHRO/0 1 *2 . *C0 S I AL PH 1 > > » 

ALST6=I ALPHO+ALPHl )/2. 

00=2.*S0RTJCHRD**2-»01**2/4. + CHRD*01*C0S « ALPHl > » 

60 TO 70 

201 WRITC(6v120>WCALC 

120 FORMATIMSHO rotor exit chokes at HAXIHUN mass flow rate = ,F9.l|> 

60 TO 2 

202 URITE(6«121) 

121 F0RMATI69M0 NO SOLUTION FOUND AFTER 100 ITERATIONS FOR CONTINUITY 
lAT ROTOR EXIT! 

60 TO 2 

203 WRITE(6»122> 

122 FORMATI60HOREOUIRED SPECIFIC WORK GREATER THAN ENERGY AVAILABLE IN 
1 GAS) 

60 TO 2 

204 WRITEI6tl23> 

123 F0RMATI74H0RE0UIRED STATOR IDEAL KINETIC ENERGY GREATER THAN ENER6 
lY AVAILABLE IN GAS) 

GO TO 2 

205 WRITEI6,124) M 

124 F0RMATJ33H0SPEC1FIC WORK REQUIRED IN SEC TOR , 1 3 , 37H GREATER THAN EN 
1ER6Y AVAILABLE IN GAS) 

GO TO 2 

206 WRITE(6tl25) H 

125 FORMAT! 55H0R0T0R IDEAL RELATIVE KINETIC ENERGY REQUIRED IN SECTORt 
113, 37H GREATER THAN ENERGY AVAILABLE IN GAS) 

GO TO 2 

207 WRITEI6,125) I 
GO TO 2 

208 WRITE16,124) 1 
60 TO 2 

209 URITE 16,128) 

128 FORMAT! 123H0THE PROGRAM CAN NOT FIND A SOLUTION SIMULTANEOUSLY SAT 
lISFYING CONTINUITY, RADIAL EO., AND THE LOSS MODEL AT THE ROTOR EX 
2IT) 

60 TO 2 
999 STOP 
END 


Blade-Row Loss Coefficient Subprograms 


The stator and rotor loss- coefficient calculations, as described in the section LOSS 
MODEL, are performed by subroutine EFFIC. Function subprograms SIMPSl, SHUB, 
and SHUB2 are used for the numerical integration required to determine the rotor hub 
wall area (eq. (C17)). SIMPSl uses a modification of Simpson's rule wherein more in- 
tervals are placed in those regions requiring greater accuracy. SHUB and SHUB2 are 
the functions being integrated. 

Program variables . - Variables transfer between main program RIFTUD and sub- 
routine EFFIC by means of both the EFFIC(ES, ER, J) arguments and labeled common 
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block /EFF/. The arguments ES and ER are the stator and rotor loss coefficients, re- 
spectively, and J is the blade- row indicator ( J = 1 for stator and J = 2 for rotor) . The 
common block variables, which were defined in the RIFTUD variable list, are ALPHO, 
ALPHl, ALUNC, BET2M, CS, CR, DO, Dl, DlA, D2M, EN, GAM, HR, HS, MU, 
NSTAR, PI, RH2RT2, SIGS, STAR, VOVCRl, W, and WOWCRM. The remaining varia- 
bles in subroutine EFFIC are defined as follows; 


APQ 

AA 

AR 

AS 

A2 

A3A2R 

A3A2S 

BB 

BWSR 

CHRDR 

D(Y) 

E(Y) 

ELOCS 

ELOSR 

ELOSS 

EWSR 

H(Y) 

K 

K2 


arithmetic statement function for Q where argument 

a 

Q 

^r 

Q 

^s 

a/2 

^3i/'^2d)^ 

(^3l/'^2D)g 

b 

arithmetic statement function for denominator of equation (10) or (11) 
where argument Y = Q or Q 

S i 

arithmetic statement function for E where argument Y = or 
(Z/c)s 

A . +A , 
w, t, r w,h,r 

arithmetic statement function for H where argument Y = or 

erix)r parameter for integration of function SHUB 
error parameter for integration of function SHUB2 


37 


I 



ref, r 


REFR 

REFS 

RER 

RES 

SRH 

SRHl 

SRH2 

SRS 

TOSR 

TOSS 

YC 



Re 


- 0 . 


9to/' 


- 0.2 


ref, s 


Re 


r 


’^r^w.h, r 

part of SRH from first integral in equation (C17) 
part of SRH from second integral in equation (C17) 

^r^ w, t, r 
(t/ s) ^ 

(t/s)g 


r 


la 


Y2 r^^-hyJ^/2 

Variables transfer among EFFIC, SIMPSl, SHUB, and SHUB2 by means of the 
function arguments and the labelled common block /SH/. The function arguments are 
SIMPS1(X1,X2, FUNC,KSIG), SHUB(X), and SHUB2(Y) where 

XI lower limit of integration 

X2 upper limit of integration 

FUNC function being integrated, SHUB or SHUB2 in this case 
KSIG error parameter, increases if integration is inaccurate 
X variable of integration 

Y variable of Integration 

The common block variables, which were defined in the EFFIC variable list, are AA, 
BB, and YC. These also are the only variables, aside from the arguments, used in 
SHUB and SHUB2. 

A further description of SIMPSl and definition of its internal variables can be ob- 
tained from reference 6. The function described is called SIMPS2, but it becomes ex- 
actly SIMPSl upon deletion of the parameter J from the function arguments. 
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Program listing , - The FORTKAN listings for subroutine EFFIC and functions 
SIMPSl, SHUB, and SHUB2 are as follows: 


SUBROUTINE EFEIC (ESfER^JI 
REAL MU 

COMHON/EFE/GANtVOVCRltWtPIfSTAR tMUfDl tALPHl « NST AR t ALPHO* ALUNC , SI6S 
1 tOOvHS^CSf UOUCRMyOI A»D2H«HRfBET2M,ENtRH2RT2«CR 
COMMON/SH/AA vBBf YC 
AfXir (GAM-1 • >/f6AH'»l. I^X^X 
D (Y >= 1 •/1.68^Y/2.88 + Y*Y/4#4^Y*Y*Y/6 .24 

HtY > = ( 1 ./X.2 + 3.^Y/1.6*5.»Y*Y/2.*7.^Y^#3/2.4^9.*Y#^4/2.8 l/OI Y> 

E (Y .92^Y/3.2*Y*V/4.8*Y^#3/6.72 >/0( Yl 
GO TO (lt2’tJ 

1 AS= AIVOVCRII 
REFS^.03734 

RES- W/PI/STAR/MU/D1/C0S(ALPH1) 

ELOCSrl. 

IFINSTAR.NE.2.0R.ALUNC.NE.0.0) E LOG S = ALUNC /2 . /S IN ( ALUNC /2 . ) 

ELOSSr ELOCS^SIGS 

A3A2S= 1 .♦(□□♦♦2-D1 + + 2 )/4./SIGS/HS/D1/EL0CS 
TOSSr .05*SI6S*STAR 

ESrCS*E ( AS >*REFS/RES*#.2«£L0SS*A3A2S/(C0S( ALPHl ) -C SA'H ( AS I ♦REFS/RES 
1«^.2*EL0SS-T0SS) 

GO TO 3 

2 ARz A(WOkCRH> 

REFR-. 11595 

CHRDRr 0.5 ♦SORT! ( 01 A -D2M--HR^HS I ( D 1A-02M ) I 

RERrw/PI*CHRDR/HR/MU/D2M/C0S( BET2M > 

ELOSR-EN*CHRDR/2./SORTI 2.I/D2M 

B«SRrPI/8.*( (D1A-02M-HR42.*HS)*(D1 A-D2M + HRI-(D1 A-02M-HR >**2) 
SRS=PI/2.^ID1A-02M-HR>*( (P I /2 • - 1 . I *0 1 A+D2M+HR ) 

EXTERNAL SHUB » SHUB 2 
K= □ 

K2r0 

A A- 4DlA-D2M-HR+2.*HSI/2. 

BB= (DlA-D2M*HR)/2. 

YC=01A/2. 

A2=A A/2. 

Y2=YC-BB/2.*S0RT(3. ) 

SRHlz2.*PI*SXMPSl(0.0tA29SHUB tK ) 

IF(K«GT.09 WRITEIGflOl K 

10 FORMAT(3H K-tI21 
SRH2-2.#PI^SlMPSl(Y2»YCtSHUB2,K2) 

IF(K 2 . 6 T .01 yRITE(6fll> K2 

11 F0RMATf4H K2-^I29 
SRHrSRHl«SRH2 
EVSR=:(SRS«SRH)/EN 
A3A2R=1.^EWSR/BUSR 

TOSR= 0.04/PI*EN^( l.-RH2RT2>/( l.^RH2RT2> 

ER=CR*EC ARI*REFR/RER^^.2*EL0SR#A3A2R/fC0S(BET2M9-CR*HIAR l♦REFR/RER 
1**.2^EL0SR-T0SR1 

3 RETURN 
END 
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FUNCTION SIMPSl (XI tX2»FUNC»KSIG) 

C THIS SOUTINE INTEGRATES FUNC(X) FROM XI TO X2 USING A MULTIPLE 

C INTERVAL SIMPSONFS RULE TECHNIQUE. 

LOGICAL SPILL 

DOUBLE PRECISION ANS.Q 

DIMENSION V (200) «H ( 200 ) « A ( 200 ) .8 (200 ) .C(200) »P (200 ) *E ( 200 ) 

DATA TWO»THREE»FOUR*THIRTy/2. 0.3.0* A. 0.30.0/ 

DATA T.NMAX.NSIG/3.0E-5. 200.1/ 

C INITIALIZE FIRST ELEMENTS OF ARRAVS. 

V(l)= XI 

H(l)= (X2 - V(l))/TWO 
A(I)= FUNC ( V ) 

B(l)= FUNC ( V(1)*H(1) ) 

C(l)= FUNC ( X2 ) 

P(l)= H(1)*(A(1) • F0UR*B(1) • C(D) 

E(l)= P(l) 

ANS= P(l) 

N=1 

FRAC=T 

SPILL=. FALSE. 

1 TEST=ABS(FRAC*ANS) 

Xi=N 

DO 3 1=1. K 

C TEST magnitude OF ATH ORDER ERROR IN THIS INTERVAL. 

IF (ABS(E(I) ) .LE.TEST) 60 TO 3 
IF (N.LT.NMAX) GO TO 2 

C GO TO FINISH IF STORAGE IS FILLED UP. 

SPILL=.TRUE. 

XSIG=KSIG*NSIG 
GO TO A 

C SUBDIVIDE INTERVAL AGAIN TO REDUCE ATH ORDER ERROR. 

2 N=N+1 

V (N)=V( I) *H( I) 

H(N) =H ( I ) /TWO 
A(N)=B(I) 

B(N)=FUNC(V(N) *H(N) ) 

C(N)=C( I) 

Pi(N)=H(N) *( A(N) ♦FOUR*B(N) *C(N) ) 

H(I)=H(N) 

B(I)=FUNC(V(I)*H(I) > 

C(I)=A(N) 

Q=P(I) 

P'( I ) =H ( I ) » ( A ( I ) .FOUR*B ( I ) *C ( I > ) 

0=P(I).P(N)-Q 
ANS=ANS*Q 
E( I )=Q 
E(N)=Q 

3 CONTINUE 

C.....TEST ALL INTERVALS AGAIN IF ANY WERE SUBDIVIDED THE LAST TIME. 

IF (N.GT.K) GO TO 1 
A- 3=0.0 

DO 5 1 = 1. N 
3=Q*E(I) 

TIGHTEN ERROR LIMIT IF TOTAL ACCUMULATED ERROR TOO LARGE. 

IF (ABS(Q/T) .LE.ABS(ANS) .OR. SPILL) GO TO b 
FRAC=FRAC/TWO 
GO TO 1 

FINISH OFF CALCULATION. 

SIMPS1=(ANS.Q/THIRTY) /THREE 
RETURN 

C THIS ENTRY USED TO GET AT INTERNAL VARIABLES. 

ENTRY SlMPXl (TT.NN.QQ) 

T = TT 

NN=N 

30=0 

RETURN 

END 
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FUNCTION SHUB(X) 

COHHON/SH/AA,BB*YC 

SHUB=(YC-BB/AA*S0RTIAA**2-X*X))*S0RTU.*BB**2/AA**2*X*X/(AA**2-X*X 
1 ) ) 

RETURN 

END 


FUNCTION SHUB2IY) 

COMHON/SH/AA,BB,YC 

SHUB2=Y*S0RT a.-»AA**2/BB**2*IY-YCJ#*2/IBB**2-tY-YCI**2) ' 

return 

END 

Rotor- Exit Continuity Subprograms 

Subroutines CONTIN and PABC provide the means for obtaining the rotor-exit con- 
tinuity solution; that is, they yield the value of ( V /V ) that satisfies the input 

\ ^ ‘^V2,m 

mass flow rate. CONTIN provides an estimate for the value of an independent varia- 
ble X that satisfies a given value of the dependent variable Y by means of curve 
fitting. The curve is a parabola whose coefficients are calculated by PABC, which is 
called by CONTIN. The estimation continues until a solution is obtained within the de- 
sired tolerance. Subroutines CONTIN and PABC are described in detail in reference 7. 

Program variables . - Variables transfer between main program RIFTUD and sub- 
routine CONTIN by means of the CONTIN (XEST, YCALC, IND, JZ, YGIV,XDEL) argu- 
ments. Transfer of variables between CONTIN and PABC is by means of the 
PABC(X, Y, A, B, C) arguments. These agruments are defined as follows-. 

XEST on call: value of X used to calculate YCALC 

on return; value of X to be used to calculate next value of YCALC 

YCALC value of Y corresponding to XEST during call 

IND on call: controls sequence of calculation in CONTIN 

on return: Indicates when a choked solution is found or when no solution can 
be found 

JZ determines whether subsonic or supersonic solution will be obtained; 

1 - subsonic solution 

2 - supersonic solution 
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YGIV value of Y desired for solution 

XDEL maximum permissable change in XEST between iterations 
X independent variable 

Y dependent variable 

A coefficient A in y = Ax^ + Bx + C 

B coefficient B in y = Ax^ + Bx + C 

C coefficient C in y = Ax^ + Bx + C 

The internal variables for CONTIN are defined in reference 7. 

Program listing . - The FORTRAN listings for subroutines CONTIN and PABC are 
as follows: 


SUBROUT I ME CONTI M( XEST ,YCALC , I ND YGIV tXaEL > 

C 

C — CONTIN CALCULATES AN ESTIMATE OF THE RELATIVE FLOW VELOCITY 
C — FOR USE IN THE VELOCITY GRADIENT EQUATION 
C 

DIMENSION XI3I«Y(3) 

NCALL = NCALL^l 

IF f IN0*NE»1 •ANO«NCALL*6T»100} GO TO 160 
GO TO < 10t30t40»S0t60«110,lS0) t INO 
C—FIRST CALL 
10 NCALL = 1 
X0RI6 = XEST 

IF t YCALC«GT.YGIV.AND.JZ»E0.1I GO TO 20 
INO r 2 
YCll = YCALC 
XIII r 0. 

XEST r XEST^XDEL 
RETURN 
20 INO - 3 

YI3I r YCALC 
XC3I = 0. 

XEST r XEST-XDEL 
RETURN 

C — SECONO CALL 
30 IND - H 

YI2I r YCALC 
XC2I = XEST-XORIG 
XEST = XEST^XDEL 
RETURN 
<10 IND r 5 

YI2I - YCALC 
XI2I r XEST-XORIG 
XEST r XEST-XDEL 
RETURN 

C — third OR LATER CALL - FIND SUBSONIC OR SUPERSONIC SOLUTION 
50 YC3I r YCALC 

XI3I = XEST-XORIG 
GO TO 70 
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60 Y 11 I r YCALC 

XIII r XEST-XORIG 

70 IF t YGIV*LT»AHlM(Yf ll,Y(2>tY<3M> GO TO I120«1301,JZ 
80 INO - 6 

CALL PABCCX,YtAPA,BPBtCPC> 

DISCR = BPB**2-<I.^APA#ICPC-Y6IV> 

IF (OISCR.LT.D. > GO TO 140 

IF CABSI4C0*^APA*ICPC-YGIV)I.LE.BPB^*2> GO TO 90 
XEST = -BPB-SIGNfSORTCOISCRI »APA » 

IF fJ2.EO*l.ANO,APA.GT.O..AND.Y(3l .GT.YCll I XEST r -BPB^ 
ISORT lOlSCRI 

IF fJZ.E0*2,AND,APA.LT,0.) XEST = -BPB-SORT ( OISCR I 
XEST r XEST/2./APA 
GO TO 100 

90 IF ( JZ.E0.2. AND.BPB.6T.0.) GO TO 130 
ACB2 = APA/BPB»ICPC-YGIV)/BPB 
IF (ABSfACB21.L£.l.E*8l ACB2rQ. 

XEST r -ICPC-Y6IV)/BPB>*I 1*^ACB2^2,^ACB2*^2) 

100 IF CXEST.GT.Xf 3M 60 TO 130 
IF CXEST.LT.Xfin GO TO 120 
XEST = XEST^XORIG 
RETURN 

C — FOURTH OR LATER CALL - NOT CHOKED 
110 IFCXEST-XORIG.GT.XC 3)) GO TO 130 
IF! XEST-XORIG.LT.Xlin GO TO 120 
YI2) - YCALC 
XC2) = XEST-XORIG 
GO TO 70 

C--THIRD OR LATER CALL - SOLUTION EXISTS* 

C--BUT RIGHT OR LEFT SHIFT REQUIRED 
120 INO = 5 
C— LEFT SHIFT 

XEST = Xf 11-XDEL*X0RI6 
XOSHFT r XEST-XORIG 
XORIG r XEST 
YI3> = YI2> 

X(31 =r XI21-X0SHFT 
YI21 r Ylll 
X<2) - Xm-XOSHFT 
RETURN 
130 IND = 4 
C--RIGHT SHIFT 

XEST = X (31+XOEL+XORlG 

XOSHFT r XEST-X0RI6 

XORIG = XEST 

Yll > - Y 121 

X<1 I = X C21-X0SHFT 

YI2I - Y<3> 

Xf2) r XI3)-X0SHFT 
RETURN 

C — THIRD OR LATER CALL - APPEARS TO BE CHOKED 
140 XEST r -BPB/2./APA 
IND = 7 

IF IXEST.LT.Xd)) 60 TO 120 
IFCXEST.GT.XI3M GO TO 130 
XEST r XEST+X0RI6 
RETURN 

C — FOURTH OR LATER CALL • PROBABLY CHOKED 
150 IF lYCALC.GE.YGIV) GO TO 110 
IND r 10 
RETURN 
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C — NO SOLUTION FOUND IN 100 ITERATIONS 
160 IND = 11 
RETURN 
END 


SUBROUTINE P ABC I X * Y , A ,B t C 1 
C 

C — PABC CALCULATES COEFFICIENTS A,B,C OF THE PARABOLA 

c — t=a*x** 2 *b*x*c, passing through The given x,y points 
c 

DIHENSION XI3>,Y(3> 

Cl - X«3)-XI1) 

C2 = IYI2»-Ylll»/«X«2I-Xill» 

A = ICl*C2-YI3l*Yll»>/Cl/iX(21-X131) 

B = C2-»XI1MX(2)1*A 

C = YU l-XIl 1*B-X(1 »**2*A 

RETURN 

END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 23, 1975, 

505-04. 
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APPENDIX A 
SYMBOLS 

A area, ft^ 

a semiaxis of ellipse in x-direction, m; ft 

b semiaxis of ellipse in y-direction, m; ft 

C loss coefficient multiplier 

Cg circumference of ellipse, m; ft 

dimensional constant, 2vr rad/ rev; 60 sec/min 
Cp dimensional constant, 1000 W/kW; 550 (ft)(lb)/(sec)(hp) 
c chord, m; ft 

Cp specific heat at constant pressure, J/(kg)(K); Btu/(lb)(®R) 

D diameter, m; ft 

E energy factor 

e blade- row loss coefficient 

g dimensional constant, 1; 32.17 (Ibm) (ft)/(lbf) (sec^) 

H form factor 

h blade or clearance height, m; ft 

Ah specific enthalpy difference, J/kg; Btu/lb 

J dimensional constant, 1; 778 (ft)(lb)/Btu 

k number of sectors 

L loss, J/kg; Btu/lb 

I surface length (leading edge to trailing edge), m; ft 

N rotative speed, rad/sec; rev/min 

Ngp specific speed, dimensionless; (ft^'^‘^) (lbm^^^/(min) (sec^^^) (Ibf^^^) 

n number of blades 

P shaft power, kW; hp 

p absolute pressure, N/m^; Ib/ft^ 

Q parameter defined by equation (12) or (13) 
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R gas constant, J/(k^(K); (ft){lbf)/(lbm)(®R) 

Re Reynolds number 

r radius, m; ft 

s blade spacing at blade- row exit, m; ft 

T absolute temperature, K; °R 

t trailing-edge thickness, m; ft 

U blade speed, m/ sec; ft/ sec 

V absolute velocity, m/sec; ft/sec 

W relative velocity, m/sec; ft/sec 

w mass flow rate, kg/sec; Ib/sec 

X X coordinate, m; ft 

y y coordinate, m; ft 

01 fluid absolute angle measured from radial direction at stations 0, 1, and la and 

from axial direction at station 2, deg 

j3 fluid relative angle measured from radial direction at stations 0, 1, and la and 

from axial direction at station 2, deg 

Y specific heat ratio 

6 displacement thickness, m; ft 

7] efficiency 

© sector angle, deg 

®cam camber angle, deg 
0 momentum thickness, m; ft 

9 

^ viscosity, (N)(sec)/m ; lb/ (ft) (sec) 

p density, kg/m^; Ib/ft^ 

$ flow angle from throughflow direction, deg 

cp stator vane stagger angle from radial direction, deg 

angle defined in figure 3, deg 
ip energy thickness, m; ft 
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Subscripts: 
av average 

b blade surface 

c clearance 

calc calculated 

cm chord mean 

cr critical condition 

df disk friction 

ex exit 

h hub 

i sector number 

id ideal 

j sector number other than mean sector 

m mean sector 

r rotor 

rad radial direction 

ref reference 

s stator 

shft shaft 

sr surface of revolution 

t tip 

tot total 

u tangential direction 

VD velocity diagram 

w wall 

X axial direction 

0 stator inlet 

1 stator exit 

la rotor inlet 


2 rotor exit 

2D two dimensional 

3D three dimensional 

Superscripts: 

* absolute total value 

” relative total value 
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APPENDIX B 


STATOR GEOMETRY MODEL 

The purpose of the stator vane geometry model is to interrelate the stator inlet and 
exit radii, inlet and exit flow angles, chord, and surface length to chord ratio in a con- 
sistent manner. The basic geometry model is shown in figure 3. There are three spe- 
cific cases considered herein; (1) an uncambered vane with known chord length, (2) a 
cambered vane with known chord length and inlet flow angle, and (3) a cambered vane 
with known inlet radius and flow angle. In all cases, the exit radius and flow angle are 
known. 

Uncambered Vane 

For this case, the vane surface length is assumed equal to the chord; that is, 

(l/c) =1. The remaining unknowns are inlet flow angle and inlet radius. From the 
s 

geometry of figure 3(a), 

r., sin <p^ 

tan <p = i ^ (Bl) 

c + rj^ cos cp^ 

Since the vane is uncambered, and = <P-^. Therefore, equation (Bl) be- 

comes 

sin 01 

tan cKq = — (B2) 

-2. + cos OL 

The inlet radius is found by using the law of cosines; that is, 

= +t\- 2cv^ cos (180 - O;^) (B3) 

Since cos = -cos (180 - a^^), we get 

rQ = + 2ct^ cos (B4) 
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Cambered Vane 


It is assumed that the surface length is a circular arc, 
arc length to chord ratio can be related to the sector angle 


as shown in figure 3 (b) . 
® . Arc length is 


The 


7 = Trr® 
180 


Chord is 


c = 2r sin — 
2 


Dividing equation (B5) by equation (B6) then yields 

1 .= 7T® 

360 sin ^ 
2 


(B5) 


(B6) 


(B7) 


The sector angle ® can be related to the camber angle ®_„^- Remembering that 
a tangent is perpendicular to the radius at the point of tangency, referring to figure 3(b), 
we can write 


(180-©^^^) +® =180 


(B8) 


Therefore, 


and equation (B7) becomes 


© = © 


cam 


(B9) 


L\ = 

c 


7T® 


cam 


© 


® 360 sin 

2 


(BIO) 


Referring to figure 3(a), the camber angle can be expressed as 


®cam = 180 > - (180 - o;^) - ^ = ^ 


(Bll) 


50 


From the law of cosines. 


cos 1 ®' = 


4 

2^1 


(B12) 


Substituting equation (B12) into equation (Bll) yields 


® = O', - Q!_ - cos 

cam 1 0 


-1 


ir^+r^-c^ 
2 ^ 0^1 


(B13) 


Equations (B13) and (BIO) are used to determine surface length to chord ratio once the 
inlet radius and chord are known. 

To determine inlet radius or chord for the two cambered vane cases being consid- 
ered, it is assumed that the stagger angle at the chord midpoint is equal to average 
flow angle; that is. 


*^cm 2 


(B14) 


Referring again to figure 3(a), the law of cosines yields 


’^cm ^ ( 2 ) ~ ^^cml o ) ^cm 


■(f)' 


(B15) 


and, since cos = -cos (180 - <iC’^^), also 


cm' 






(B16) 


Adding equations (B15) and (B16) gives 




(B17) 
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Known chord. - In this case, we are trying to find the inlet radius. With the 
chord known, equation (B15) is solved for r„_ using the quadratic formula 

OlLL 


r =i 
cm 2 


c cos + 
^ ctn 


it 


c cos (p 




(B18) 


With r^jjj evaluated thusly, equation (B17) yields the inlet radius 



(B19) 


Known inlet radius . - In this case, we are trying to find the chord. Subtracting 
equation (Bl^ from equation (B16) and solving for r^^ result in 


r 

cm 



2c cos cp 


cm 


(B20) 


Substituting equation (B20) into equation (B17) yields 



(B21) 


Using the quadratic formula and taking the positive root finally give 


c = 


(B22) 
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APPENDIX C 


ROTOR GEOMETRY MODEL 

The purpose of the rotor geometry model is to provide a consistent basis for eval- 
uating mean surface length, chord length, blade-surface and end-wall surface areas, 
and number of blades. Shown in figure 4 is the rotor geometry with dimensions ex- 
pressed in terms of various radii and the stator blade height. The model is based on 
the meridional-plane projection of the rotor. The tip contour is assumed to be circular 
and the hub contour is assumed to be elliptical, each being a 90° arc. It is further as- 
sumed that the rotor consists of full blades only. 


Surface and Chord Lengths 


The mean surface length and mean chord are the elliptical arc and straight line, 
respectively, joining the mi(%)oints (see fig. 4) of the rotor leading and trailing edges. 
The general formula for the circumference of an ellipse with semiaxes a and b is 


C 


e 



(Cl) 


Therefore, the mean surface length of the rotor is 




r 




From the Pythagorean theorem, the chord length is 



(C2) 


(C3) 
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Wall Areas 


The blade surface area is the area between the elliptical hub and the circular tip. 
The general formula for the area within an ellipse is 

A = Trab (C4) 

Therefore, the blade surface area for one passage (i. e. , two sides) is 




’^2, t V <^la - ^2, h> - (^la ’ *'2, t/ 


(C5) 


The tip and hub wall areas are determined as the surfaces of revolution of the tip 
and hub curves around the tuibine axis. The general formula for the surface of revolu- 
tion of a curve y = f(x) between (x^^, and (x 2 , is 



With the turbine axis taken as y = 0, the centers of the tip and hub curves are at x = 0 
and y = Therefore, the equation of the circular tip is 


+(y - 




(C7) 


while that of the elliptical hub is 


2 

a 



(C8) 


where the semiaxes a and b are 


a 


^la - ^2 


t"'^s 


(C9) 


and 
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(CIO) 




For the circular tip, differentiation of equation (C7) yields 


dx y - ’^la 


(Cll) 


Substitution of equations (Cll) and (C7) into equation (C6) results in 


^sr,t = 27T(ria- ^2,t> 


^’^la-’^2,t 


1 

1 

1 


dx 


(C12) 


Integrating equation (C12) and dividing by the number of rotor passages then gives the 
tip wall area for one passage 


V,t, r ^ (^la“’^2,t^ 


--l)^la^^2.t 


(CIS) 


For the elliptical hub, the area of the surface of revolution cannot be analytically 
expressed because equation (C6) cannot be integrated analytically . Further, neither 
form of equation (C6) can be integrated numerically over its entire range because of an 
infinite slope at one limit. However, if we break the curve into two pieces, we can 
write 



From equation (C8) , the two derivatives can be determined as 
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bx 


(C15) 





and 


^ bVb2-(y-rig)2 

Substituting equations (C8), (CIS), and (C16) into equation (C14) and dividing by tbe 
number of rotor passages then give the following hub wall area for one passage: 



Equation (C17) can be integrated numerically. In the program, each of the integrals is 
evaluated by Simpson's rule. 


Blade Number 


From the approach of reference 8 along with the additional assumption that 
V 1 a ~ U- a number of blades can be computed as 

U, ±ci 

n^ = 2 tt tan (C18) 
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This number of blades is based on not permitting the velocity to fall below zero any- 
where within the rotor channel; the limiting condition, therefore, is zero velocity on 
the pressure surface at the rotor inlet. Equation (C18) yields a rather large number of 
blades, especially at the higher flow angles. Studies such as the one reported in ref- 
erence 5 have shown that the number of blades can be reduced significantly from these 
high values without significant degradation in performance. Therefore, a reduction 
factor is used in equation (C18). Also, the small difference between and is 

neglected. With these changes, the number of rotor blades is computed herein as 

n^ = -2- (110 - a )tan oi^ (C19) 

30 ^ ^ 
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TABLE I. - INPUT FORM WITH SAMPLE DATA 


STATEMtNT 

number 


FORTRAN STATEMENT 


. THIS IS A SAMPLE CASE USING S , I UWITS 

$ IN PUT R = 20 8.1, GAM= 1 , 6 67 , ALPHA 0 = 5 5 . 6,, A L PH A 1 = 7 2,. 0,.N = 4 
PTIN=9.101,W=,277 1 , P0W = 22 . 371, RIRIA =1.0216, R0R1A=1 .2 
RH2RT2 =0_. 34 9 3, O PT 2 = . 0 023 , M U= . 58 04 E- 0 4 , K= 5 , ru=l , NSTAr 
SHIS^ THE^ USING U.S. ^ C USTO MARY UNIT 

$ INPUT R=3 8 . 68 , N= 38500 . , TTIN=19 5 0 . , PTIN = 13 . 2 , W=[. 611 , 

1^4 4 4- . 4- 4--4-F-f— H — I — F-’-h ) 1- +-H — — 4_ 4 — i -4_4. _4_4 — 4 -? 1 — | — 4 — 4—4 - +Z_4-. 4_ _| — , — 4 —, — ^ 


p 51 , 

5 81,^ RT 

= 0 I 


TTI,N^,ip 
2R1A=^0 . 


83., ,5,5^ 

715.5.. . 


P0W=5 0^, ,MU= ,^5 9E- 04 , ^ 




4- f. 4- _-4-- » 4 - 4- -+ — ( 


-1 1-- 1 ^— + ^ -f +- -I— f 4- V 4- ) 


t - ♦ 1—4)44 


- 4 - ) » 4-4 -4 4--F- 


)- ) 4- -4 — 4 — 4 -F-H 


I 2 3 U 5 I 6 I 7 8 9 10 II 12 1 15 It 15 16 l7 18 1 19 20 21 22 25 2t|25 26 27 28 29 30| 31 52 35 3t 35 36|37 38 39 tO U t2|t3 UUUS U6 A7 A0|t9 50 51 52 53 5<4|55 56 57 58 59 60|6l 62 63 6t 65 66|67 68 69 70 71 72| 73 74 75 76 77 78 79 8 

NASA-C-836 (REV. 9-14-59) 






TABLE n. - OUTPUT FOR SAMPLE CASE USING SI UNITS 


RADIAL INFLOW TURBINE VELOCITY DIAGRAM DESIGN ANALYSIS 
THIS IS A SAMPLE CASE USING S5 UNITS 

THIS OUTPUT IS IN THE FOLLOWING UNITS. 


TEMPERATURE 

KELVIN 


PRESSURE 
N/SO CM 


GAS CONST 
JLS/KG-K 


ROT SPEED 
RAD/SEC 


MASS FLOW 
KG/SEC 


POWER 

KW 


VISCOSITY 
N-SEC/SO M 


VELOCITY 

H/SEC 


SPEC WORK 
JLS/6N 


DIAMETER 

CM 


ANGLE 

DEGREES 


♦INPUT* 



INLET TEMP = 

1083.3300 

STATOR IN angle = 55.60 

INLET PRESS= 

9.1011 

STATOR EX angle = 72.00 

ROTAT SPEEOr 

9031.710 

STATOR ASPECT RATr .0000 

MASS FLOW = 

.2771 

OIAM RATIOS 

SHAFT POWER= 

22.371 

STAT IN/ROT IN =1.2581 
STAT EX/ROT IN =1.0216 
ROT EX TP/ROT IN= .7133 
ROT EX HUB/TIP = .3993 


GAS CONSTANT r208.U10 
SPEC HT RATIOrX.6670 
VISCOSITY r .580M-0«I 


♦OUTPUT* 


STATOR KE LOS COEFrl.OOQO 
ROTOR KE LOSS COEFrl.OOQO 
ROTOR IN/OEL RVU =1.0000 
ROTOR EX RAD SECTS= 5 

ROT EX SECT/HN RVUrl.OO 1.00 l.OQ 1.00 1.00 


CL HT/ROT EX TIP D= .0023 




ABS 

ABS 



ABS 

ABS 

ABS 


REL 

REL 

REL 

REL 

REL 


DIA- 

TOTAL 

TOTAL 

STATIC 

STATIC 

flow 

VELO- 

CRIT 

BLADE 

TOTAL 

TOTAL 

FLOW 

VELO- 

CRIT 


METER 

TEMP 

PRESS 

TEMP 

PRESS 

ANGLE 

CITY 

VEL RAT 

SPEED 

TEMP 

PRESS 

ANGLE 

CITY VEL RAT 

STATOR INLET 

19.550 

1083.33 

9.101 

1072.02 

8.865 

55.60 

108.99 

• 209 







LOSS COEFr , 

0690 








stator 

HGT= 1.8527 

.NUMBER 

OF VANES 

= 16.0 


STATOR EXIT 

15.875 

1083.33 

8.995 

1013.99 

7.623 

72.00 

268*67 

• 506 







ROTOR INLET 

15.539 

1083.33 

8.995 

1010.85 

7.565 

71.92 

279.59 

.517 

313.29 

1020.95 

7.746 

-31.50 

99.93 

• 194 

LOSS COEF= . 

2820 














ROTOR EXIT 

3.872 

926.12 

5.925 

908 .32 

5.699 

• 00 

136.09 

.277 

78.05 

931.98 

6.019 

-29.83 

156.86 

.319 


9.593 

926.12 

5.919 

908.99 

5.699 

.00 

133.99 

.272 

92.59 

939.36 

6.096 

-39.74 

162.96 

.330 


6.035 

926.12 

5.886 

910.66 

5.694 

• 00 

126.79 

• 256 

121.66 

990.35 

6.115 

-93.83 

175.69 

.355 


7.978 

926.12 

5.852 

912.83 

5.699 

• 00 

117.60 

.290 

150.79 

997 .96 

6.203 

-52.09 

191.19 

• 385 


8.920 

926.12 

5.810 

915.93 

5.699 

• 00 

105.95 

.215 

179.82 

957.21 

6.310 

-59.61 

208.96 

• 918 


10.363 

926.12 

5.762 

918.50 

5.699 

• 00 

69.06 

.181 

206.90 

968.07 

6.936 

-66.91 

227.09 

.953 


11.089 

926.12 

5.735 

920.20 

5.699 

• 00 

78.98 

• 160 

223.99 

979.11 

6.507 

-70.65 

236.82 

.970 



MASS 

OlAG 













DIA- 

FLOW 

SPEC 

TOTAL 

STATIC 











METER 

RATE 

WORK 

EFFIC 

EFFIC 

NUMBER 

OF ROTOR 

BLADE$= 

12.0 







3.872 

.0000 

81.769 

.919 

.839 











9.593 

.0919 

81 .769 

.916 

• 839 











6.035 

• 0516 

81.769 

.907 

• 839 











7.978 

.0592 

81.769 

• 896 

• 839 











8.920 

.0631 

81.769 

.883 

.839 











10.363 

• 0617 

81 .769 

• 868 

• 839 











11.089 

• 0000 

81.769 

• 860 

.839 











♦OVERALL performance* 

TOT-TOT PRESSURE RATIO = 1.5553 
TOT-STAT PRESSURE RATIO= 1.6125 
DIA6 AV6 SPECIFIC WORK = 81.7686 
DJA6 TOTAL EFFICIENCY = .8915 
DIAG STATIC EFFICIENCY = .8339 
NET TOTAL EFFICIENCY r .8800 
NET STATIC EFFICIENCY = .8232 


LOSS/IDEAL T-S DEL M 
STATOR = .0252 
ROTOR = .0779 
WINDAGE = .0098 
CLEARANCEr .0059 
EXIT KE = .0658 

SPECIFIC SPEEDr .739 


TABLE in. - OUTPUT FOR SAMPLE CASE USING U. S. CUSTOMARY UNITS 


fHlS IS THE SAME CASE USING U.S. CUSTOMARY UNITS 


THIS OUTPUT IS IN THE FOLLOmING UNITS, 


TEMPERATURE 

PRESSURE 

GAS CONST ROT 

SPEED 

MASS FLOW 

POWER 

VISCOSITY 

VELOCITY SPEC 

WORK 

diameter 

angle 

DEG RANKINE 

LB/SO IN 

BTU/LB-R 

REV/MIN 

LB/SEC 

HP 

lb/ft-sec 

ft/sec BTU/LB 

IN 

DEGREES 

♦INPUT* 















INLET TEMP = 

1950.0000 

STATOR IN ANGLE : 

55.60 

STATOR 

KE LOS 

COEFrl.OOaO i 

GAS CONSTANT 

: 38.6800 



INLET press: 

13.2000 

STATOR EX ANGLE : 

72.00 

ROTOR 

KE LOSS 

COEFrl .0000 

SPEC HT RATI0:1.6670 




ROTAT speed: 

38500.000 

STATOR ASPECT RAT: 

.0000 

ROTOR 

IN/DEL RVU :l.DDDO 

VISCOSITY 

: .3900-04 



-ass flow : 

.6110 

DIAM 

RATIOS 



ROTOR 

EX RAD SECTS: 5 







SHAFT power: 

30.000 

STAT IN/ROT 

IN :l 

.2581 

ROT EX 

SECT/MN 

RVU:l .00 

1.00 1 

.00 1.00 1.00 






STAT EX/ROT 

IN :l 

• 0216 












ROT 

EX TP/ROT in: 

.7133 












ROT 

Ex HUB/TIP : 

.3493 

CL HT/ROT EX TIP D= .0023 






♦OUTPUT* 

















ABS 

ABS 



ABS 

ABS 

ABS 


REL 

REL 

REL 

REl 

PEL 


DIA* 

TOTAL 

TOTAL 

STATIC 

STATIC 

FLOW 

VELO- 

CRIT 

BLADE 

TOTAL 

total 

FLOW 

VELO- 

CRIT 


METER 

TEMP 

PRESS 

TEMP 

PRESS 

ANGLE 

CITY 

VEL RAT 

SPEED 

TEMP 

press 

ANGLE 

CITY 

VEL RAT 

STATOR INLET 

7.697 

1950.00 

13.200 

1929.63 

12.858 

55.60 

355.93 

.204 







LOSS COEf: . 

0640 








STATOR 

HGT: .7294 

, NUMBER 

OF VANES 

: 16.0 


STATOR EXIT 

6.250 

1950.00 

13.046 

1825.10 

11.057 

72.00 

881.46 

.506 







ROTOR INLET 

6.118 

1950.00 

13.046 

1819.53 

10.973 

71.92 

900.88 

.517 

1027.71 

1836.81 

11.235 

-31.50 

327.85 

. 19 4 

LOSS COEF: . 

2820 














ROTOR EXIT 

1.524 

1667.02 

8,593 

1634.98 

8.186 

.00 

446.49 

.277 

256.36 

1677.56 

8.729 

-29.83 

514.70 

.319 


1.808 

1667.02 

8.577 

1636.19 

8.186 

• 00 

437.98 

.272 

303.76 

1681.86 

8.769 

-34.74 

533.00 

.330 


2.376 

1667.02 

8.537 

1639.23 

8.186 

.00 

415.83 

.258 

399.16 

1692.64 

8.869 

-43.83 

576. 4C 

.355 


2.944 

1667.02 

8.487 

1643.09 

8.186 

.00 

385.83 

.240 

494 .56 

1706.34 

8.996 

-52.04 

627.26 

• 385 


3.512 

1667.02 

8.427 

1647,78 

8.186 

• 00 

345.96 

.215 

589.96 

1722.97 

9.152 

-59,61 

683.91 

.418 


4.080 

1667.02 

8.357 

1653.30 

8.186 

• 00 

292.19 

.181 

685.36 

1742,53 

9.335 

-66.91 

745. C5 

.453 


4.364 

1667.02 

8.318 

1656.37 

8.186 

• 00 

257.47 

.160 

733.06 

1753.41 

9.438 

-70.65 

776.96 

.470 



MASS 

OIAG 













DIA- 

FLOW 

SPEC 

TOTAL 

STATIC 











METER 

RATE 

WORK 

EFFIC 

EFFIC 

NUMBER 

OF ROTOR blades: 

12.0 







1.524 

.0000 

35.178 

.919 

.634 











1.608 

• 0914 

35.178 

.916 

.834 











2.376 

.1138 

35.178 

.907 

.834 











2.944 

• 1305 

35,178 

• 896 

.834 











3.512 

.1392 

35.178 

.883 

.634 











4.080 

• 1361 

35.178 

.868 

.834 











4.364 

• 0000 

35.178 

• 860 

.834 











•OVERALL PERFORMANCE* 

TOT-IOT PRESSURE RATIO = 1.5553 
TOT-STAT PRESSURE RATI0= 1.6125 
01A6 AVG SPECIFIC UORK = 35.1777 
0IA6 TOTAL EFFICIENCY = .8915 
D1A6 STATIC EFFICIENCY r .8339 
NET TOTAL EFFICIENCY = .8800 
NET STATIC EFFICIENCY = .8232 


LOSS/IDEAL T-S DEL H 
STATOR = .0252 
ROTOR = .0779 
WINDAGE = .00M8 
CLEARANCEr .0059 
EXIT KE : .0658 

SPECIFIC SPEED: 95.355 


O 











